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Neutrinos, and primarily neutrino oscillations, have undoubtedly been one of the most 
exciting topics in the field of high-energy physics over the past few years. The existence 
of neutrino oscillations would require an extension of the currently accepted description 
of sub- nuclear phenomena beyond the Standard Model. Compelling evidence of new 
physics, which seems to be pointing towards neutrino oscillations, is coming from the 
solar neutrino deficit and from the atmospheric neutrino anomaly. More controversial 
effects have been observed with artificially produced neutrinos. The present experimental 
status of neutrino oscillations is reviewed, as well as the planned future experimental 
programme, which, it is hoped, will solve most of the outstanding puzzles. 

1. Introduction 

Since the discovery of the electron neutrino in 19550, the study of neutrino physics 
has often led to great progress in the unravelling of sub-nuclear phenomena. How- 
ever, although they have been used as powerful probes of elementary particles in- 
teractions, very little is known about the general properties of neutrinos: we know 
neither their mass, nor their magnetic moment; we don't know whether they are 
Dirac or Majorana particles; we don't know how many neutrino species exist in 
nature. In many respects, neutrinos perhaps remain the most intriguing mystery in 
particle physics. 

In more recent years most of the research activity in neutrino physics has fo- 
cused on the phenomenon of neutrino oscillations, originally postulated by Bruno 
PontecorvoH Interest in this field has been revived by several very interesting experi- 
mental results, which in some cases very strongly support the hypothesis of neutrino 
oscillations. This is only possible if, contrary to what assumed in the Standard 
Model of electro-weak interactions, neutrinos have non-zero and non-degenerate 
masses. Thus, evidence for neutrino oscillations would represent a major discovery 
in elementary particle physics and would then require a complete revolution in our 
understanding of the natural world. 

This is a very exciting time for neutrino physics, with new data becoming avail- 
able in a relatively short period of time. Whilst this paper was being completed, 
the long-awaited first results from the SNO experiment were presented.B We hope 
that many others will follow, from all the running or planned experiments, and that 
more data will soon lead to the solution of this fascinating puzzle of particle physics. 
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In this paper, after a brief introduction to massive neutrinos and to the formahsm 
of neutrino osciUations, the latest experimental results will be reviewed. The focus 
will be on those experiments which have added pieces of information to the current 
picture of neutrino oscillation phenomenology, either by observing a signal or by 
setting a limit in an interesting portion of the relevant parameter space. Planned 
experiments, designed to resolve some of the as yet unsolved ambiguities, will also 
be described. 

2. Neutrino Masses 

The present experimental limits on neutrino masses, as obtained from direct 
kinematic searches, under the assumptixjiLjJiat flavour and mass eigenstates are 
essentially coincident, are the following D'EHstD 

m^^ < 2.2 eV (95% C.L., from ^^Hc + + 9^)] 

m^^ < 170 keV (90% C.L., from tt+ ^ n+ + v^) (1) 



< 15.5 MeV (95% C.L., from r ^ 57r(7r°) + i/^) . 



Although the smallness of these limits might lead to the naive conclusion that 
neutrino masses are exactly vanishing, there is no theoretical reason to assume that 
this is indeed the case. On the contrary, it is generally believed that particle masses 
can be identically equal to zero only if they are associated with an exact gauge 
symmetry. This is not true for the lepton quantum number, which is then expected 
not to be conserved. 

Perhaps the most popular expression for the neutrino mass term is the one 
associated with the so-called see-savS form of the mass matrix. In the simple case 
of one generation this mass term is given by: 



-) Z) [Z) • 

where a right-handed neutrino singlet, vr, has been introduced. Assuming that the 
Dirac mass mo is much smaller than the Majorana mass Mm? the diagonalisation 
of (|^) gives a heavy eigenstate, of mass ^ Mm, and a light eigenstate, which can 
be identified with the standard left-handed neutrino and whose mass is given by: 

/ON 

If the Dirac mass is of the same order as that of a fundamental fermion, be 
it a quark or a charged lepton (rriD ~ 0(raq,i)), and if Mm ~ 0{Mgut), then 
the neutrino mass turns out to be naturally very small, and in a range of values 
compatible with those indicated by experimental results. 



"in the generic n— generation case m^) and Mjv/ are n X n matrices: then m^j <C A/^j means that 
all eigenvalues of mj5 are small compared to all eigenvalues of Mjvf . 
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Several versions of the see-saw mechanism are considered in literature.cl Among 
them, the so-called quadratic see-saw, in which all the Majorana masses are similar 
to each other, and therefore the light neutrino masses scale as the Dirac masses 
squared: 

Another possibility is the linear see-saw, where Mm oc mo and then: 

: : ~ iT^qi,ii ■ "^q^j^ ■ rriq^^g^ , (5) 

Both versions imply a hierarchy between the light neutrino mass eigenstates: 

m^^ <C TO^^ <C m^^ . (6) 



3. Neutrino Oscillations in Vacuum 

Ordinary neutrinos are produced in weak processes together with their charged 
partner leptons (W~^ — > i'^i'i), as for example in pion decays or in nuclear beta 
decays, and they are therefore eigenstates of the weak interactions. On the other 
hand, weak interaction eigenstates and mass eigenstates are not necessarily coin- 
cident. With a change of basis, neutrinos produced as flavour eigenstates can be 
represented as a coherent linear superposition of the mass eigenstates: 

l^a) k.) (a = e,M,T) (7) 

i 

where i runs over all the existing mass eigenstates and Uai are the elements of a 
unitary matrix.* 

Although from measurements of the width at LEP we know that the number 
of light active neutrinos is Afu — 2.994±0.012l3, nothing forbids the existence of more 
than three neutrinos of definite mass. If, for example, there were four independent 
mass eigenstates, a change of basis from the mass to the flavour eigenstates would 
necessarily lead to one linear combination, 

i^,) = ^c/,: k,), (8) 

i 

which has no normal weak coupling. The linear combination is therefore referred 
to as a "sterile neutrino" , as opposed to the standard "active" neutrino. 

As we shall see in the next sections, the existence of at least one sterile neutrino 
might be the only way to accommodate all the experimental results on neutrino 
oscillations. Nonetheless, for the sake of simplicity, in the following discussion 
we will only consider the three neutrino flavours whose existence has been firmly 

^ The convention of using Greeks letters to identify flavour eigenstates and Latin letters to identify 
mass eigenstatates has been adopted throughout this paper. 
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established.'^ However, it must be noted that the same formahsm would be valid for 
any number of neutrinos. 

After a neutrino has been prepared in the state | i^q ) at i ~ 0, each of its 
components will evolve according to the Schrodinger equation. At a time t after 
production it will be: 

I ^„ (i > 0)) = J2 u:, ) = E E ^/^^ e"''''* ^» I ) ' (9) 

i 13 i 

Thus, after the neutrino produced in the flavour eigenstate a has travelled a 
distance L t, the probability of finding it in a different flavour state /? is non-zero 
and is given by: 

- ^p) - l(^MOko>P = lE^Up^ e-^"* c/;/ . (10) 

Using the unitarity condition for U, after some algebra the probability (|l^) takes 
the following oscillatory form: 



Vil^a Vp) 



Am%L 



where Amf^ — \m1 — m?j \ and the ultra-relativistic expansion for the neutrino energy 
has been used.'' Thus the transition probability V{va — > vp) depends on a combi- 
nation of the mixing matrix elements, on the squared difference between the mass 
eigenvalues of the mass eigenstates, and on two parameters, the neutrino energy E 
and the distance L between the production and the observation points, which are 
determined by the experimental conditions. 

Although from the simplicity of the preceding discussion it might seem that 
neutrino oscillations can be obtained as a straightforward consequence of elementary 
quantum mechanics, a rigorous treatment of the problem would require that a wave 
packet formalism be considered. We shall not dismss this approach here and the 
interested reader is remanded, for example, to ref.O and refs. therein. 

The lepton mixing matrix U , also known as the Maki-Nakagawa-Sakata (MNS) 
niatrix,E3 is the leptonic analogue of the CKM matrix for the quark sector. In the 
case of three neutrino flavours the explicit relation between the flavour and the mass 
eigenstates is given by: 



(12) 




'^the existence of the r neutrino, of which only indirect evidence ™as available until very recent 
times, has been proven by the emulsion-based DONUT experiment, t3 at Fermilab. 

Ei = Jp^ + ni1 ~ p + ml /2p ~ p + ml /2E. 
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It is sometinies convenient to parametrize U in a form anafogous to that of the 
CKM matrixjj which for Dirac neutrinos is: 

C12C13 S12C13 sise""'^ 

U = \ -S12C23 - ci2S23Si3e*'^ C12C23 - si2S23Si3e*'^ S23C13 | , (13) 

S12S23 - Cl2C23Sl3e** -C12S23 - Sl2C23Sl3e*^ C23C13 

with the standard convention aj — cosOij and Sij — sinOij — 1,2,3 are the 
generation labels). Thus in this formalism the mixing matrix depends only on three 
mixing angles, 9i2, 9i3 and ^23, and a CP-violating phase, (5/ 

Although the most general description of the mixing among all neutrinos implies 
the three-flavour treatment, most experimental results are expressed in the two 
generation mixing representation. More recent analyses, however, generally attempt 
a global fit to all data assuming three-flavour mixing. 

In the case of transitions between two flavours the mixing matrix has the form 
of a rotation in a two-dimensional space: 

... / cos6' sin0\ 

^=[-sine cos^j' (14) 

where 9 is the mixing angle, analogous to the Cabibbo angle for the quarks. The 
relation between the flavour and the mass eigenstate is then given by 



\iya) ~ COS 9\vi) + sm 9\V2) 
\vp) = - sin 1 1^1 ) + cos 1 1^2 ) ' 

The oscillation probability in this case can be written as follows: 

L 



(15) 



V{va -> i^p) = sin^ 29 sir? (tt- — ) (16) 

where Xosc is the oscillation wavelength, defined as: 

AttE 

Xosc = -r — 2 ' ^^"^f 

where E is again the neutrino energy and Am^ is the squared neutrino mass differ- 
ence. 

Substituting the expression of Xosc in (p^, the oscillation probability takes its 
usual form: 

Vii^c. - lyp) = sir? 29 sin' j^ ^-^^Am^^ ^ ^^^^ 

where Am? is expressed in eV^, the detector-to-source distance L is in meters (A:m) 
and the neutrino energy E in MeV (GeV). 



^in the case of Majorana neutrinos there would be two additional phases. 
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The oscillation probability (18) is the product of two factors. The first one, 
sin^ 29, does not depend on the experimental conditions, but is an intrinsic parame- 
ter which describes the strength of the coupling between the two neutrino flavour; it 
gives the maximum amplitude of the neutrino oscillations. The mixing is maximal 
when sin^ 29 = 1, namely when 9 — tt/A. 

The second factor is another oscillatory term, the period of which is determined 
by Am^ L/E. In order to have sensitivity to small values of Am?, the ratio E/L 
must be also small. This can be achieved by having low energy neutrino beams and 
large distances between the neutrino source and the detector. On the other hand, 
when the oscillation phase becomes too large, that is to say when E/L is too small 
compared to Arri^, the oscillations occur very rapidly and, due to the finite energy 
resolution of the experiments, the term sin^(l. 27 Am^ L/i?) averages to 1/2. 

Before concluding this overview of the phenomenology of neutrino oscillations 
in vacuum, it is interesting to consider the effect of the CPT, CP and T symmetries 
on the oscillation probabilities. Without entering in any details, we shall only say 
that for both Dirac and Majorana neutrinos, CPT invariance implies: 

Vii^c f/s) ^ V{D(3 Va) (CPT conserved). (19) 

In particular, this means that the neutrino and the antineutrino survival prob- 
abilities are the same: 

Viva^ Va) ^'Pi^a^ Va) (CPT conscrvcd). (20) 

On the other hand, the oscillation probabilities for neutrinos and antineutrinos 
are generally different. They coincide only if CP, the symmetry which converts a 
left-handed neutrino into a right-handed antineutrino, is conserved: 

V{va i^is) = 'Pi'^a i^p) (CP conserved). (21) 

Finally, if the time-reversal symmetry T is conserved, the oscillation probability 
is invariant under interchange of the initial and the final states: 

'P(t^a -> J//3) = Vivfj -^Va) (T conserved). (22) 



4. Neutrino Oscillations in Matter 

In the previous section we have seen that the probability for neutrino oscillations 
in vacuum cannot exceed sin^ 29, 9 being the mixing angle in vacuum between 
the two flavours under consideration. On the contrary, the situation can be very 
different for neutrino oscillations in matter. Matter effects can greatly enhance 
neutrino mixing, resulting in a large oscillation probability even for very small 
vacuum mixing angles. 

When they travel through matter, neutrinos of all flavours can have neutral- 
current interactions with the protons, neutrons and electrons of the medium. How- 
ever only electron neutrinos can interact with the electrons, undergoing a coherent 
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forward scattering via a W boson exchange. The consequence of this asymmetry 
hetweeii neutrino flavours is known as the Mikheyev-Smirnov-Wolfenstein (MSW) 
ejfecMB 

At low neutrino energies, for electron, muon and tau neutrinos traversing an elec- 
trically neutral and unpolarised medium, the matter- induced potentials are given 
by: 



2 



(23) 

2 



where Gf is the Fermi constant and Me and are the electron and the neutron 
number densities respectively. For antineutrinos the potentials have opposite signs. 

For simplicity, the discussion will be restricted to the two-flavour scenario, more 
specifically to the mixing between and v^. fn this case, using the ultra-relativistic 
expansion for the neutrino energy and neglecting all common phases (which would 
not affect the result for mixing between active neutrinos), the time evolution for 
the flavour eigenstates is given by: 



d 
i— 
dt 



I 



V 



COS 26* + ^GpMe sin 20 



— — sm 20 COS 20 

AE AE 



(24) 



where the common term due to neutral-current interactions, proportional to the 
identity matrix, has been ignored. 

The simplest case of matter distribution is that of constant matter density. Me = 
const. Although not very realistic, this is however quite an instructive example and 
we shall discuss it here in some detail. 

As for the mass and flavour eigenstates for oscillations in vacuum, the mass and 
flavour eigenstates in matter are also connected by a two-dimensional rotation. The 
relation between the two bases is given by: 

) = cos Om ) + sin Om \ v^, ) , , 

= -sin^Mke) + cosOm\v,). ^ ' 

The vacuum mixing angle, 0, has been replaced by the mixing angle in matter, Om, 
given by: 



sir 

\ 2E J 

■ cos 20 - ^f2GFMe) + sin^ 261 



sin^ 261 



sin^ 2Qm ^-^^ 5 • (26) 

2E 



ft has to be noticed that, since the mixing angle is not the same as in vacuum, the 
matter eigenstates | ) and | v^^ ) are not coincident with the mass eigenstates in 
vacuum, I V\ ) and ^2 ) • 
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The oscillation probability is of the same form as in Eq.(p^), but the vacuum 
mixing angle and the vacuum oscillation wavelength are now replaced by those in 
matter: 



Vii^e J^m) = sin' sin' ( tt^ ) , (27) 



M 



with Xm given by: 



2tt , , 

Am = , , , = . (28) 

cos 20 - ^2G,M^ \{^\" sm^ 2B 



\ 2E I \ 2E 

From Eq.(^^ it follows that, regardless of the smallness of the mixing angle 
in vacuum, the mixing angle in matter can be very large. In particular, maximal 
mixing can be achieved if the medium density is such that the following condition 
is satisfied: 

\/2G'fA/; = ^^cos26'. (29) 
2E 

This is the the so-called MSW resonance condition. Using the characteristic solar 
electron density. Me ^ 10^^ cm~^, and the characteristic value for solar neutrino 
energies, E ^ 1 MeV, assuming small mixing in vacuum (hence cos 26' ~ 1), one 
obtains Am^ ^ 10^^ eV^. 

Since A/'e > 0, Eq.(|29|) is fulfilled only if Am^cosfl > 0. Once a convention on 
the phase has been chosen, Eq.(|2^) implies that resonant oscillation enhancement 
is possible only for one particular sign of Ato^. For example, if cos 6* > 0, neutrino 
oscillations are enhanced if Am^ — m\ — m\ > 0. On the other hand, from the 
change of sign in the matter-induced potentials, it follows that for antineutrinos the 
resonance condition would require Am^ < 0. This means that, for a given sign of 
Am^, matter effects cannot enhance neutrino and antineutrino oscillations at the 
same time: if neutrino oscillations are enhanced, antineutrino oscillations will be 
suppressed, and vice-versa. 

The evolution equation ( p4[ ) cannot be solved analytically for any non-uniform 
matter distribution. In general a numerical solution of the system of differential 
equations ( p^ ) has to be computed, accounting for the density profile and the neu- 
trino energy distribution. An interesting case to consider is the so-called adiabatic 
approximation for matter density monotonically decreasing along the neutrino path. 
Without entering into the details of the calculation, we shall only qualitatively dis- 
cuss the results for the case of two-flavour neutrino mixing (for example between v^. 
and f^), which is useful to describe neutrino oscillations in the Sun. 

Electron neutrinos produced in the core of the Sun, where the density is well 
above that corresponding to the MSW resonance, will initially see a mixing angle 
Ocore ~ tt/2. From Eq.(^) it follows that, at production point, a neutrino born as a 
Ve will essentially coincide with one of the matter eigenstates 1^2 ^ . The adiabaticity 
condition, that requires slowly changing matter distribution, guarantees that the 
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neutrino system can gradually adjust to the changing density of the environment 
and therefore the system will not make any transition to the other matter eigenstate. 
As 1^2^ propagates through the mantle, it encounters regions of smaller densities: 
the effective mixing angle decreases and the strength of the mixing increases, until 
the resonance condition is fulfilled and maximal mixing is reached. As the neutrino 
travels further, the mixing angle becomes smaller and smaller, approaching the 
value of the mixing angle in vacuum, 9. If 6* is very small, the Ve component of 
is very small at the final point and 1^2^ is mainly composed of t^^. The survival and 
the oscillation probabilities for electron neutrinos are given by: 



-> i^e) = sin^ ( 



(30) 



Thus, in the range of validity of the adiabatic approximation and for small values 
of the mixing angle in vacuum, if the depth of the traversed matter is large enough, 
the probability of finding the neutrino in the flavour state Vf, when it gets outside 
the Sun is tiny and a complete conversion of v,, to u^^ is possible. 

At this point it becomes natural to discuss the solar neutrino deficit and to 
present the status of the search for neutrino oscillations in the Sun. 

5. The Solar Neutrino Problem 



5.1. Solar Models 

Historically the first hint for neutrino oscillations came from the observation of 
neutrinos from the Sun. In 1968 Ray Davis andJiis collaborators published the first 
results of the Homestake chlorine experiment, Fir^ showing that the measured flux 
of solar electron neutrinos was significantly lower than the expected value. This 
was the beginning of the solar neutrino deficit saga. 

The Homestake results_ss£ere later confirmed by founpther experiments: SAGEjl3 
GALLEX,S Kamiokandejli and Super-Kamiokande@ Both SAGE and GALLEX 
are radiochemical experiments, like Homestake, while Kamiokande and Super-Kamio- 
kande are water-Cherenkov experiments. Before we discuss them, it is important 
to introduce the so-called Standard Solar Model (SSM), used to predict the solar 
neutrino fiux under the assumption that no exotic phenomena, such as neutrino 
oscillations, can affect the nature of the produced neutrinos before they reach the 
Earth. 

Our Sun belongs to the category of the main sequence stars, which produce 
energy in their interiors via thermonuclear reactions. The two reaction chains re- 
sponsible for the energy production are known as the pp-cycle and the CNO-cycle. 
For both cycles the overall result is the fusion of hydrogen nuclei into helium, with 
the emission of electron neutrinos: 



Ap+2e 



^ He -I- 2z^e + 26.73 MeV. 



(31) 
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Only ^ 2% of the solar energy is emitted in the form of neutrinos, while the rest is 
radiated through photons. 

In the three decades since the first experimental results on solar neutrinos, there 
has been a flourishing of solar models, all of which use the present values of some 
fundamental solar parameters (such as radius, mass, luminosity, He/H ratio) as 
constraints on the stellar evolution. They all make essentially the same assumptions 
on the nature of the solar energy and of the solar energy transfer mechanism. The 
most famous and accredited model is certainly that of Bahcall-Pinsonneault (BP 
SSM)xiJ In the following we shall be referring to its two most recent versions, the 
BEgsEland BP200o|l in terms of which most published experimental results are 
interpreted. All reactions belonging to the pp and the CNO cycles, together with 
the corresponding average and maximum neutrino energies are listed in Tab.|l| (from 
ref.E). The CNO cycle contributes less than 2% to the total neutrino flux. 

Table 1. Average and maximump4iputrino energy for the different reactions contributing to the 
solar neutrino flux (tab. from refo). 



Source Reaction < > (MeV) Max. E^ (MeV) 



PP 


P + 


P 


d + e+ - 


-f Ve 


0.2668 


0.423 ± 0.03 


pep 


P + 


e" 


" +p ^ d- 


+ 


1.445 


1.445 


hep 


3He + 


P 


^ 4H<= + e 


+ +Ve 


9.628 


18.778 


7 Be 


e~ - 




Be Li 


+ 


0.3855 


0.3855 












0.8631 


0.8631 


8B 


8B- 




*Be* + e+ 


+ 


6.735 ±0.036 


~ 15 


13n 


13n 




>i3c + e+ 


+ Ve 


0.7063 


1.1982 ± 0.0003 


15 o 


150 




>i^N + e+ 


+ V£ 


0.9964 


1.7317 ±0.0005 


17-p 






.1^0 + e+ 


+ V£ 


0.9977 


1.7364 ±0.0005 



The neutrino fluxes predicted by BP2000, as well as the predicted neutrino 
capture rates in the chlorine and gallium experiments, are listed in Tab.|. The 
uncertainties on the neutrino absorption cross sections represent the dominant con- 
tribution to the error on the gallium rate predictions, while for the chlorine rate 
those uncertainties are relatively small compared to other systematic components 
intrinsic to the experimental technique. The rates are measured in Solar Neutrino 
Units (SNU), a convenient unit to describe the rates of solar neutrino experiments 
(1 SNU = 10-36events atoms-^s-^). 

The energy spectm^of neutrinos from the pp-cycle are shown in Fig.^ (BP2000 
predictions, from refnj). The uncertainties on the different flux components are 
also shown and the energy thresholds of the different experiments are indicated at 
the top of the plot. 

There are several reasons, other than the high precision of the input data and 
calculations, why the BP SSM is believed to be robust. First of all, there is what 
is known as the "luminosity constraint" argument. It assumes that the Sun is in 
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Neutrino Energy (MeV) 



Fig. 1. BP2000E3 predictions for the energy spectra of solar neutrinos produced in the pp-cw^Je 
reactions, which are responsible for more than 98% of the energy balance in the Sun (fig. fron£3). 
The continuum spectra are expressed in events cm^'^s^^MeV"^ at one astronomical unit, while 
the monochromatic lines (pep and ^Be) are given in events cm^^s"'^. 



0.1 
0.08 
0.06 
0.04 
0.02 


-0.02 
-0.04 
-0.06 - 
-0.08 
-0.1 



Sound velocities 



Bahcall-Pinsonneault 98 
LOWLl + BISON Measurement 



♦-'Be lowered {Icr off Ga) 



Fig. 2. Comparison between the predicted, (Model) and measured (Sun) sound speeds in the Sun. 
Predictions come from thaJ3P98 modelJ£3 while measurements come from the helioseismological 
data on the sound speeds .PH 
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Table 2. BP200oF^ predictions for the solar neutrino fluxes and the neutrino capture rates in 
chlorine and gallium experiments. 



Source Flux (lO"' cm^^g-i) CI (SNU) Ga (SNU) 



PP 



13 



5.95(l.OOlo - 69.7 

pep 1.40 X 10-2(l.00lg:°i5) 0.22 2.8 

hep 9.3 X IQ-'^ 0.04 0.1 

''Be 4.77 X 10-1(1. OOt°;}o) 1.15 34.2 

®B 5.05 X 10-*(l.OOto ig) 5.76 12.1 

N 5.48 X 10-2(1. ooto;?7) 0.09 3.4 

l^O 4.80 X 10-2(1. OOto ig) 0.33 5.5 

i^'F 5.63 X 10-*(l.00t|3J^) - 0.1 

Total 7.6+J-5 128+^ 



a stationary state and then infers a strong correlation between the solar neutrino 
flux and the solar luminosity (£© = 2.4 x 10"^^ MeV s-^, known to an accuracy of 
^ 0.4%). From Eq.(pl|) we see that, in the nuclear reactions occurring in the Sun, 
the production of two neutrinos is accompanied by an energy release of 26.73 MeV. 
Neglecting the small amount of energy carried away by the neutrinos, this means 
that the electron neutrino flux in the absence of oscillations can be estimated as 

where Rq is one astronomical unit (= 1.496 x 10^^ cm). 

Another model-independent constraint, although looser than the luminosity one, 
is given by the fact that the helium nuclei, necessary for the boron and beryllium 
production, are created in the pp and pep reactions. Thus: 

(^Be) + («B) < (pp) + (pep) . (33) 

One the most impressive checks of the BP SSM come perhaps from the com- 
parison between the BP98 model predictions, obtained without any adjustment of 
the PEimmeters, and the most accurate helioseismology measurements of the sound 
peed.E^ Fig.|| shows the excellent agreement between the calculated and the mea- 
sured values: the size of the fractional difference between measurements and predic- 
tions (0.001 standard deviation for radii between 5% and 95% of the Sun radius) is 
much smaller than any generic change in the model with an impact on the predicted 
neutrino flux. 



S" 
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5.2. Experimental Results 

We shall now describe the experimental techniques which have been used to de- 
tect solar neutrinos. The experiments can be subdivided into two categories: rad io- 
chemical experiments (Homestake, SAGE, GALLEX and its successor, the GNOcJ 
project) and water-Cherenkov experiments (Kamiokande, Super-Kamiokande and 
SNO). One of the main differences between these two categories is the fact that, 
while for radiochemical experiments only an integrated measurement of the fluxes is 
possible, water-Cherenkov detectors can perform real-time measurements and there- 
fore can study differential distributions (i.e. energy spectrum, angular distribution, 
correlation with th e S un position in the sky). 

The HomestakelialiS chlorine experiment is located 1480 m underground (4200 
m.w.e., meters water equivalent), at the Homestake gold mine. Lead, South Dakota, 
USA. A tank of 6 x 10^ litres volume is filled with tetrachlorethylene (C2CI4). The 
detection principle is based on the reaction: 

i/e +^^G\ ^e- + ^^Ar(Ti/2 = 35 d) , (34) 

which has a threshold of 814 keV, above the end-point of the pp-cycle energy spec- 
trum. After an exposure time of 1 — 3 months, the '^^Ar atoms produced in ( p^ 
are extracted by purging the detector with ^He and detecting the Auger electron 
produced in the electron capture of the radioactive '^^Ar nuclei. The measured 
rate, which mainly arises from of andjBe neutrinos, plus a small contribution 
from pep neutrinos, is 2.56 ± 0.23 SNUjij significantly below the predicted value 
of 7.6t.i'\ SNU. Although the Homestake detector is not calibrated, for lack of a 
suitable artificial neutrino source, the efficiency of the extraction technique has been 
checked by doping the detector with a known amount of radioactive argon atoms. 

All gaUium experiments, SAGE,E1 GALLExH and GNO,izl detect neutrinos 
from the Sun by means of the process: 

Ve Ga ^ e- Ge(Ti/2 = 11.43d) , (35) 

for which the energy threshold (233.2 keV) is well below the maximum energy of 
the pp neutrinos. This means that gallium experiments, by measuring low-energy 
solar neutrinos, have the ability to prove the thermonuclear nature of the energy 
production mechanism in the Sun. 

The GALLEX experiment, now discontinued as such, was located at the Gran 
Sasso Underground Laboratory (LNGS), Assergi, Italy. It consisted of 30.3 tons of 
^^Ga, in the form of a concentrated solution of gallium chloride (GaCla — HCl) in 
water. The ^""^Ge atoms form the volatile compound GeCU which, at the end of 
each run (3 — 4 weeks) is swept out of the solution by means of a nitrogen stream. 
The GeCU is absorbed in water in a gas scrubber, where the nitrogen is filtered out, 
and then converted to GeH4. This compound is finally introduced, together with 
xenon, into a proportional counter where the number of ^^Ga atoms is determined 
by counting the radioactive decays. The final results from GALLEX(I-IV)Ea give 
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a solar neutrino flux of 77.5^7 g SNU, again well below theoretical predictions of 
1281? SNU. The GALLEX detector was calibrated using two independent methods. 
The first one made use of neutrinos from two intense (> 60 PBq) ^^Cr sources: the 
combined value of the ratio R between the neutrino source strength as derived from 
the measured rate of ^^Ga production and the directly determined source strength 
is i? = 0.93 ± 0.08.Ea The other calibration was performed by injecting, under vary- 
ing conditions, a known amount of ''^As into the full-scale detector .lj The arsenic 
isotope decays by electron capture and positron emission to ^^Ge {T1/2 — 2.72 d), 
producing radioactive atoms which mimic the solar neutrino capture kinematics. 
Although neutrinos from the ^^Cr source provide a better match to solar neutrino 
energies, the second method has the great advantage of large statistics. The mea- 
sured recovery rate of ^^Ge from gallium is 1 ± 0.01. The two results@@ together 
rule out the possibility that the solar neutrino deficit observed by GALLEX can be 
attributed to systematic errors in the radiochemical extraction procedure. 

Th. SAgS experiment is e„ned o„. by a R„s.i.„-A„„ica. »llabo,...o. 
at the underground Baksan Neutrino Observatory (4715 m.w.e.), in the North- 
ern Caucasus mountains. The detector, which weighs 57 tons, uses gallium in its 
metallic form. The germanium produced in (^5|) is removed from the metallic gal- 
lium by a liquid-liquid extraction into a HCl — H2O2 phase. After this initial step, 
the experimental technique is very similar to that followed by GALLEX. The re- 
sults reported by the SAGE collaboration are based on a 10 years exposure time 
(from January 1990 to Dctober 1999), yielding a measured solar neutrino flux of 
74l7;4(stat.-|-sys.) SNU,ll3 in very good agreement with the GALLEX results. The 
SAGE detector has been calibrated, exposing 13 tons of metallic gallium to an in- 
tense 51 Cr source.E_l As for GALLEX, the ratio of the measured production rate to 
that predicted from the source activity, R = 0.95 ± 0.11(stat.)lQ'Q8, has confirmed 
that the extraction procedure is reliable. 

In 1998, after an upgrade of the existing experimental setup, the GNO (Gallium 
Neutrino Observatory) experiment0 started operation at LNGS, as the successor 
project to the GALLEX experiment. The main aims of GNO are to provide a long 
time record of low-energy solar neutrinos, to determine the bulk production rate 
with an accuracy of 5 SNU and to monitor the time dependence of the pp-neutrino 
flux during a whole solar cycle to a precision of ^ 15%. The GNO collaboration 
has recently reported the results from the GNO I phase of the experiment, based 
on 19 months of observation (from May 199E.^to January 2000). The measured solar 
neutrino flux is 65.811q 2(stat.-|-sys.) SNUj23 which confirms the GALLEX result. 
A combined analysis of GNO I and GALLEX gives 74.11^;^ SNU.0 

Both KamiokandeB and Super-Kamiokande@ are water-Cherenkov detectors, 
located 1000 m underground (2700 m.w.e.), in the Kamioka mine in Japan. 

Kamiokande^started its operation in 1984 and was originally intended to study 
nucleon decay.Ed It was later upgraded to detect also low energy events and suc- 
ceeded in observing the first solar neutrinos in 1987. Kamiokande was made of a 
cylindrical tank, of total volume 4.5 kton, filled with pure water. An inner volume of 
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2.14 kton was defined by 980 inward-looking photo- multiplier tubes (PMTs), used 
to detect the Cherenkov light produced by the relativistic particles traversing the 
water. Kamiokande ended physics data-taking at the beginning of 1995 and was 
completely stopped in summer 1997. Since then its physics programme has been 
continued by its successor experiment, Super-Kamiokande. 

Super-Kamiokande consists of a huge cylindrical tank (volume 50 kton) filled 
with pure water (water transparency is ~ 100 m at A = 420 nm). Also in this case 
an inner volume (16.9 m diameter, 36.2 m height, volume 32.5 kton) is defined by 
an inner surface equipped with a large number of PMTs (11, 146 PMTs of 50 cm 
diameter, 40% surface coverage). The increased coverage of PMTs makes Super- 
Kamiokande not only bigger than Kamiokande, but also better in terms of energy, 
position and angular resolution. An outer volume, 2 m thick, equipped with 1, 185 
PMTs (20 cm diameter) surrounds the inner detector and serves as an active veto 
counter against gamma-rays, neutrons and through-going cosmic muons. For the 
solar neutrino analysis, the fiducial mass is 22.5 kton, the fiducial volume boundaries 
being 2 m inside the inner surface. 

In both Kamiokande and Super-Kamiokande the solar neutrinos are detected 
through the observation of the Cherenkov rings produced by the electrons emitted 
in the elastic process: 

+ ^ i^x + e~ (x = e, /i, r) . (36) 

Although the cross-sections for this process are very small, ^elastic interactions turn 
out to be very useful, thanks to the nice correlation between the recoil electron 
momentum direction and the direction of the incoming neutrino [Ef,Orec < 2me, 
6 ~ 18" for 10 MeV neutrinos). Therefore the direction of the recoil electron can 
be used to correlate the direction of the impinging neutrino with the Sun's position 
in the sky (see Fig.^) . Moreover the energy of the recoiling electron can be used to 
obtain a lower limit on the incoming neutrino energy. 

The total released energy is correlated to the number of PMT hits, which is in 
turn a function of the total emitted light, corrected for light absorption through 
water, PMT geometry and overlapping hits. The energy scale, however, must be 
determined in an independent manner. 

Kamiokande used gamma-rays from the Ni(n,7)Ni reaction to obtain it.tj The 
measurement precision was at the 1 — 2% level, the main limitation being the knowl- 
edge of the branching ratios and of the neutron absorption cross-sections for different 
nickel isotopes. Only limited information could be extracted about the energy res- 
olution, while no information at all could be obtained on the angular resolution of 
the detector. 

Super-Kamiokande determines the neutrino energy scale by injecting, at various 
positions of the inner detectors, electrons from an electron LINAC placed near the 

f 0.920 X 10~*''(' — — — ) cm^ for electron neutrinos and 0.157 X 10~^^( — — — 1 cm^ for muon 
VlOMeV/ VlOMeV/ 

neutrinos. 
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Fig. 3. Super-Kamiokande (1258 days) distribution of the cosine etithe angle 9sun between the re- 
coil electron momentum and the vector from the Sun to the EarthE3 The dots are the experimental 
data, while the solid histogram is the best fit for signal plus background. The peak emerging at 
cos 8sun over a flat background is due to solar neutrinos. 



detector tank.E3 Electrons from the LINAC span the energy range 5—16 MeV, 
which matches exactly the region relevant for solar neutrinos. By means of this 
procedure, the absolute energy scale is known to better than 1%. Moreover, Monte 
Carlo studies show that the energy resolution of the detector can be reproduced to 
an accuracy of 2%, while the angular resolution is reproduced to better than 1.5° 
for 10 MeV electrons. 

In order to check these results and, what is possibly even more important, in 
order to reduce the intrinsic systematic error inherent in the described procedure, f 
another calibration technique has been introducedBil A pulsed deuterium-tritium 
neutron generator produces isotropically distributed neutrons of 14.2 MeV, via the 
reaction '^H H He + n. The neutrons create ^^N by the {n,p) reaction on 
-^^O in the water. ^^N decays (t = 7.13 s, Q — 10.4 MeV) are characterised by 
the emission of a 4.3 MeV electron in coincidence with a 6.1 MeV gamma-ray and 
they are therefore well suited to verify the solar neutrino absolute energy scale. 
Additional ^^N is also naturally produced in cosmic muon capture by ^^O. These 
events, which are selected out of the full data sample and which are generally used 
to monitor the solar neutrino signal extraction method, can be also analysed to 
check the absolute energy scale. Excellent agreement (within 1%) has been found 
between the data and the MC, previously tuned on the LINAC data. Moreover 
it has been possible to show that the position and the angular dependence of the 
energy scale are within systematic uncertainties for the energy scale of the detector. 

The energy threshold for Kamiokande and Super-Kamiokande is determined by 

f electrons from the LINAC only move in a downward direction; due to limited access, the beam 
pipe and the LINAC calibration equipment can only be operated at a finite number of positions; 
the presence of the beam pipe in the tank limits the calibration procedure precision and is in fact 
the main source of systematic uncertainties at low energy. 
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the threshold for detecting the recoil electron in the elastic scattering (pq), which 
is about 5 — 7.5 MeV. Thus the solar neutrino flux measured in those experiments 
is essentially the *B flux (plus a small component from the hep reaction in Super- 
Kamiokande). The value measured by Kamiokande, based on 2079 days of data- 
taking (from January 1987 to February 1995), is 2.80 ± 0.19(stat.) ± 0.33(sys.) x 
10^ cm~^s"^.E3 The Super-Kamiokande collaboration has recently reported a result 
based on 1258 days of data-taking, for recoil electrons in the energy range 5 — 
20 MeV. The measured neutrino flux is 2.32±0.03(stat.);° °^(sys.) x 10^ cm-^s-ifl 
For the hep neutrinos they set an upper limit of 4 x 10** cm~^s~^, at 90% C.L. 

Before we move to other issues relevant for the understanding solar neutrino 
puzzle, we summarize the results discussed up to now and their implications on 
neutrino physics. In Tab.^ the measured values of the_integrated solar neutrino 
fluxes are compared with the BP2000 model predictionaH^. 

Table 3. Comparisoip-between the measured integrated solar neutrino fluxes and the corresponding 
BP2000 predictionsO Units are SNU for chlorine and gallium experiments and lO^cm -^s ^ for 
and hep neutrinos. 



Experiment Measured BP2000 Measured/BP2000 



Chlorinfli 


2.56 ± 0.23 




0.34 X (l.00±0.06) 


GALLEX+GNoBi] 


74 1+''-^ 


128t? 


0.58 X {l-OOto lD 


SAGE0 


75.4l^:» 


1281? 


0.59 X (1.00 ±0.12) 


* B-Kamiokande0 


2.80 X (1.00 ±0.14) 


5.05 X (l.OOlo ig) 


0.55 X (l.00l°it) 


* B-Super-KamiokandeEi 


2.32 X (1.00 ±0.03) 


5.05 X (l.00l°;?°) 


0.46 X (l.OOlo;?g) 


/iep-Super-Kamiokande0 


< 40 X 10-3 


9.3 X 10-3 


< 4.3 



In brief, all experiments see a deficit of solar neutrinos compared to predictions. 
We have seen that the solar model's predictions appear to be robust, and so are the 
experimental results (although the Homestake experiment is sometimes still being 
questioned because of the lack of a calibration). Therefore neither an astrophysical 
solution to the solar neutrino puzzle, nor attributing the discrepancy between the- 
ory and experiments to not completely understood systematic uncertainties seem 
acceptable. There must be a particle physics solution to the solar neutrino deficit, 
which is however more than a sole normalisation problem. The impossibility of 
accommodating all the experimental results in a consistent picture forces us to go 
beyond the standard physics framework and to look for explanations which involve 
new phenomenology. 

The internal inconsistency of the solar neutrino results is very nicely summarised 



18 An Experimentalist's View of 



Total Rates: Standard Model vs. Experiment 

Bahcall-Pinsonneault 98 




Theory ■ 'Be □ p-P- pep Experiments b 

»B ■ CNO 



Fig. 4. Measured vs. predicted solar neutrino rates (fig. from ref.E3) for all the solar neutrino 
experiments. The predicted rates for the various components of the neutrino spectrum are taken 
from the BP98 modehB 



in Figjj (from ref , where a breakdown of all the measured and calculated event 
rates is shown as a function of the experimental technique and of the different re- 
actions contributing to the solar neutrino flux. First of all there is the low rate 
observed in the chlorine experiment, which basically measures '^B neutrinos, with 
smaller contributions from ^Be and CNO neutrinos. Second, there is the incompati- 
bility between the chlorine and the water-Cherenkov experiments. In fact, since the 
shape of the ^B neutrino energy spectrum is very stable under reasonable changes 
of the solar parameters, in the absence of new physics it should be possible to 
compute the ^B rate in CI from the spectrum observed by Kamiokande and Super- 
Kamiokande, which essentially measure the ^B flux at high energies. But the num- 
ber obtained from the analysis of the water-Cherenkov data is by itself higher than 
the total flux measured by Homestake. This result, taken as such, implies that the 
net contribution of the pep, ^Be and CNO reactions to the total chlorine experi- 
ment flux is negative. In other words, the Homestake measurement is compatible 
with the ^Be neutrinos being completely suppressed, and this is indeed a big puzzle, 
since beryllium is needed for the production of ^B neutrinos. Finally, examining the 
gallium experiments result in detail, one finds that essentially all of the measured 
flux can be accounted for by only considering the contribution from the pp-cycle, 
which is known within 1% and therefore is not questionable. Ignoring for simplicity 
the small contribution from ^B neutrinos (which, by the way, can also be derived 
from water-Cherenkov results), this means that there is no room to accommodate 
the significant contribution expected from the ^Be lines. 
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The strong correlation between the predicted ^Be and the fluxes is shown in 
Fig.^ for many different solar models, together with the best fit solution for the ''Be 
flux as obtained from experimental data. All the fluxes have been normalised to the 
predictions of the BP98 model and a Sql, rectangle has been drawn, showing that 
all models but one (that of Dar-ShavivO) agree with each other within 3cr of the 
BP98 model and are all far away from the value derived from experimental data. 
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Fig. 5. Correlation between the '''Be and tiw *B fluxes (normalised to BP98 predictions) as given 
by published solar models (fig. from ref.tj). The rectangle corresponds to the 3a contour for 
the BP98 model prediction and the ellipses are the 1, 2, 3cr contours for the best fit solution as 
obtained from all experimental data assuming no new physics. 

In summary, if one believes the conclusions of the standard solar models, the 
solar neutrino puzzle can only be solved if one is either ready to discard at least 
three experimental results (the chlorine experiment and either the gallium or the 
water-Cherenkov measurements), which seems unreasonable, or to accept that new 
physics, not accounted for in the standard model of electro- weak interactions, is 
modifying the energy spectrum of solar neutrinos. 

One very attractive possibility is that neutrino oscillations are changing the 
neutrino flavour composition while they travel between the centre of the Sun and 
the Earth. If solar electron neutrinos have transformed into muon or tau neutrinos 
before reaching the detector, given their low average energy they would fail to 
produce muons or taus via charged-current interactions. A "disappearance" of 
electron neutrinos in the chlorine and gallium experiments would then become an 
expected effect. The deficit in the water-Cherenkov experiments, where in principle 
muon and tau neutrinos could be detected via their neutral-current scatters with 
electrons of the medium, could also be explained by the low cross-section for that 
process. 

Interpreting the solar neutrino deficit in the two-fiavour mixing framework, al- 
lowed regions can be drawn on the oscillation parameter space {sin^ 26, Am^} at a 



20 An Experimentalist's View of 




Fig. 6. Allowed regions (99 % C.L.) on the plane (sin^ 2f?, Am^), as obtained in rcf. from 
a combined two-flavour analysis of the total rate results in the chlorine, SAGE, GALLEX and 
Super-Kamiokande experiments, for oscillations in matter (left) and in vacuum (right). 

set C.L..'' For example, Fig|| shows the 99 % C.L. solutions as obtained in rcfS 
from a combined analysis of the total rate results in the chlorine, SAGE, GALLEX 
and Super-Kamiokande experiments, for oscillations in matter (MSW solutions, 
left) and in vacuum ("Just-so" or VO solutions, right). The matter solutions are 
conventionally named according to the allowed ranges of the oscillation parameters: 
small mixing angle (SMA, Am^ 10"^ eV^, sin^ 26* 10"'' - 10"^), large mix- 
ing angle (LMA, Am^ - 10"^ - 10"'' eV^, sin^ 26* > 0.5) and low Am^ (LOW, 
Am^ ~ 10"'^ eV^, sin^ 29 > 0.9). The VO solution (Am^ < 10"^ eV^) corresponds 
to the case when one astronomical unit exactly matches the neutrino oscillation 
length for typical solar energies. 

One way to discriminate between small and large mixing angles is to consider 
the so-called "day/night effect". Neutrinos detected during night, as opposed to 
neutrinos detected during day, are subject to matter effects in the Earth, which are 
expected to be sizeable for large values of the mixing angle. The solar neutrino flux 
can be studied as a function of the solar zenith angle 9z and an asymmetry A can 
be defined: 

where $jv f^nd are the night and day fluxes respectively. A non-zero asymmetry 

'^As stressed mref.i , under certain conditions the MSW solutions to the solar neutrino deficit can 
extend to the tt/4 < 8 < tt/2 domain, which is not covered by the standard two-flavour analyses; to 
overcome this problem, results can be more conveniently presented in the {tan-^ 0, Am^} parameter 
space. 
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would favour large mixing angles. Another way to separate between the two regimes 
is to consider the shape of the solar zenith angle distribution, which is expected to 
be flat for large mixing angles, while for small mixing angles it should show an 
excess of events in the upward-going direction (around — 1 < cos < —0.8). 
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Fig. 7. Super-Kamiokande (1258 da*«|) measured distributions for the solar zenith angle (left) and 
of the solar neutrino energy (right). e3 
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Fig. 8. Super-Kamiokande (1258 days) measurement of the seasonal variatiptt of the solar neutrino 
flux, not corrected for the effect due to the orbit eccentricity of the Earth.E3 



Oscillations would also manifest themselves as a distortion of the neutrino energy 
spectrum and it is therefore important to study the neutrino energy distribution, 
which for Super-Kamiokande is approximated by the recoil electron energy. It turns 
out that, if the predicted hep neutrino flux is correct, a distribution peaked at high 
energies would favour the VO solution. 

Finally, strong evidence for neutrino oscillations in vacuum would come from 
the observation of a seasonal variation of the solar neutrino flux, different from the 
standard geometrical one due to the orbit eccentricity of the Earth(e = 0.0167).* 

' For Am^ of the order of that suggested by the VO solution, and for solar neutrino energies of the 
order of ~ few MeV, the seasonal change in the distance between the Sun and the Earth (AHq = 
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The day-night asymmetry obtained by SupetJCamiokande based on 1258 days 
of data is A = 0.033 ± 0.022(stat.)lo:oi2(sys.),Ei which is only a 1.3cr deviation 
from zero asymmetry. The zenith angle distribution (Fig.^ left) is rather flat and 
shows no feature which might lead to a preference for the SMA solution. The 
measured energy spectrum shows no statistically significant distortion relative to 
the predicted ®B spectrum (Fig.^, right). No deviation from the expected standard 
annual variation of the neutrino flux has been observed (see FigJ|) . 

The Super-Kamiokande collaboration has recently performed a flux-independent 
analysis of their data, assuming two-flavour mixing. The resulting exclusion plot 
(95% C.L.) for the Ve v^^r case is shown in Fig.|| (shaded areas), together with the 
allowed areas obtained using the zenith angle distribution and the SSM predictions 
(dotted) and the allowed areas obtained from a combined analysis of the GALLEX, 
SAGE, Homestake and Super-Kamiokande flux measurements (hatched). From this 
analysis the SMA and the VO solutions seem to be disfavoured at 95% C.L. The 
same analysis repeated for pure v^, — > Vs shows that all solutions are disfavoured at 
95% C.L. for the sterile neutrino hypothesis. 



-3 




-4 -3 -2 -1 

10 10 10 10 , 1 

sin^(20) 



Fig. 9. Superr+jamiokande exclusion areas for solar neutrino oscillations (shaded area, based 
on 1258 days).E3 The allowed areas obtained using the zenith angle distribution and the SSM 
predictions (dotted line) and those obtained combining all the experimental results on integrated 
solar neutrino fluxes (hatched islands) are also shown. 

Although Super-Kamiokande has measured the energy dependence and the zenith 
angle dependence of the solar *B flux, and extracted useful information about the 

2ei?Q) is a sizeable fraction of the oscillation length in vacuum, A = iirE/l^m? (numerically, 
2 X 0.016 X (1.496 x 10" m) x (2 x lO"!" eV^)/l MeV ~ 1). 
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MSW and the VO solutions out of an impressively high-quality and high-statistics 
sample, the solar neutrino problem can by no means be considered definitively 
solved. Many questions still remain open and it is therefore fair to say that still 
lot of information is missing and further experimental effort is req uire d. Along 
this line there come the Sudbury Neutrino Observatory project (SN0),c3 which has 
been running for about one year, as well as the BorexinoEj and the KamLANDcj 
experiments, which are both scheduled to start data taking very soon. 

SNO is a 1 kton heavy-water Cherenkov detector located 2 km underground 
(6010 m.w.e), in the Creighton mine, Sudbury, Ontario, Canada. The detector is 
made of a spherical acrylic vessel, 12 m diameter, containing ultra-pure D2O and 
surrounded by an ultra-pure H2O shield, in turn contained in a cylindrical cavity 
(34 m height, 22 m maximum diameter). The light produced in the water is detected 
by 9,456 PMTs (20 cm diameter) with light concentrator, installed on a stainless 
steel structure surrounding the acrylic vessel. 

As well as using the elastic scattering (ES) reaction (|3^), as Kamiokande and 
Super-Kamiokande do, SNO can detect solar neutrinos via charged-current (CC) 
and neutral-current (NC) interactions on deuterium: 

Ve + d^p + p + e- (threshold = 1.4 MeV) (38) 

and: 

Vx + d^n + p + Vrc (threshold = 2.2 MeV) . (39) 
For completeness, we give again the ES reaction (^6|): 

i^x + J^x + e~ {x = e, \i, t) . (36) 

Neutral current interactions do not depend on the flavour of the incident neutrino 
and therefore neutrino oscillations to active neutrinos should not affect the number 
of NC events. On the other hand, given the energy range for solar neutrinos, and 

coming from oscillated Vf, would be below threshold for the production of their 
charged partners: the CC sample would then be depleted by neutrino oscillations. 
This means that the CC/NC ratio can be used as a powerful probe to test the 
oscillation hypothesis. Moreover, following the same argument, any distortion in 
the measured CC energy spectrum with respect to predictions would be another 
signature of new physics. 

The electron produced in (^), as opposed to that produced in (|36|), carries 
most of the neutrino kinetic energy in the final state. Therefore it provides a better 
estimate of the initial neutrino energy (resolution of ^ 20% for the range of interest) 
and can be used to observe possible distortions of the neutrino energy spectrum. 

While electrons produced in CC and ES reactions are detected via the asso- 
ciated Cherenkov light, neutrons from reaction ( ^9|) are detected using different 
techniques during the three different phases of the experiment. In the first phase, 
recently concluded, when only heavy water was used, they were measured through 
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the 6.25 MeV photons produced in the neutron capture (25% efficiency); in the 
current phase, when 2.5 tons of NaCl have been added to the water, they are seen 
through the 8.6 MeV photons produced in the neutron capture on CI (85% effi- 
ciency); finaUy, and this wiU be SNO's third phase, the sah will be removed and 
^He proportional counters will be installed, allowing direct detection of the neu- 
trons (45% efficiency) and measurement of the spectra with completely different 
systematics. 

The sample collected during the first phase, when the sensitivity to NC was 
lower, is essentially dominated by CC interactions, with a smaller component from 
ES events. Therefore the CC/NC comparison will not be possible until data from 
the second phase, with the NC rate highly enhanced by the addition of salt, be- 
come available. However some information on neutrino oscillations can already be 
extracted from the first phase data, by comparing the CC and ES rates. In fact, 
while the CC reaction ( ^8|) is only sensitive to i^e's, because of the low solar neu- 
trino energy, the ES reaction ( ^ ) is sensitive to all active neutrinos, although with 
a reduced sensitivity to and Vr- Thus, if neutrino oscillations change the flavour 
composition of the solar neutrino flux, transforming a fraction of the original i^e's 
into other active neutrino species, the CC sample will be depleted, while the ES 
sample remains unchanged. In other words, if the CC and ES measurements of the 
®B electron neutrino fluxes give: 

$fB^(z.e)<<I>FB^(^e), (40) 

then, without having to refer to any particular solar model, the inequality between 
the two measured values can by itself be interpreted as an indication of neutrino 
oscillations to active neutrinos. 

The SNO collaboration has very recently published their first results on ES and 
CC reactions, based on 240.95 live days of data.H In order to avoid biases in the 
event selection procedure, the data has been sub-divided in two subsets, with 70% 
of them used to establish the analysis procedure and the remaining 30% used as 
a blind sample to validate the analysis itself. The final data set consists of 1169 
events. 

The different contributions from CC, ES and residual neutron events are ex- 
tracted by means of a maximum likelihood function, which combines the informa- 
tion from three variables: the effective kinetic energy of the event, the angle between 
the reconstructed direction of the event and the instantaneous Sun-to-Earth direc- 
tion and a volume-weighted radial variable, related to the position of the event 
inside the detector. The breakdown of the events as obtained from the likelihood 
maximisation is given in Tab.^. 

SNO data show that the shape of the CC energy spectrum is consistent with 
the Bahcall SSM predictions within one standard deviation. However, the low ratio 
of the measured flux to SSM predictionsBj 0.347 ± 0.029J3 confirms the solar 
neutrino deficit. Moreover it is smaller than what obtained by Super-Kamiokande 
from the analysis of ES events.H 
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Table 4. Contribution from CC, ES and neutron events to the neutrino interaction sample ex- 
tracted from 240.95 live days of SNO data, first phasej^ The quoted errors are statistical only. 



Sample 


Nr. of events 


CC 


975.4 ±39.7 


NC 


106.1 ± 15.2 


neutron 


87.5 ± 24.7 



Using events with kinetic energy above 6.75 MeV, the *B fluxes as measured 
from CC and ES interactions in SNO are: 

'^SNoi'^e) = {1.75 ± Om{stat.)tnll{sys.) ± 0.05{theor.)) x 10^ cm^^g-i 

(41) 

^SNoi^x) = (2.39 ± OM{stat.)to[ll{sys.)) x 10^ cm^^s-i (42) 

where the theoretical uncertainty in (41) is that on the total CC cross section. 

The difference between ^sNoi^e) and '^sNoi^^) is (0.64±0.40) x 10^ cm^^g-i^ 
that is to say a discrepancy of only 1.6 <j. However the SNO measurement of the ES 
flux is consistent with that measured with higher precision by Super-Kamiokande:L3 

$f|(i/^) = (2.32 ± 0.03(sta<.)t°:°^(st/s.)) x 10^ cm'^s'^ , (43) 

which permits to make a direct comparison between ^g^o('^<^) ''^^'^ ^sxi'^x)- 

The difference between the ES flux measured by Super-Kamiokande and the 
CC flux measured by SNO is (0.57 ± 0.17) x 10^ cm~^s~^, this time corresponding 
to a more significant 3.3 a discrepancy. This result, by showing that the elec- 
tron neutrino flux from CC interactions is indeed smaller than that obtained from 
ES events, provides quite a convincing hint for appearance of non-electron active 
neutrino flavours in the solar neutrino flux. 

If electron neutrinos oscillated to sterile neutrinos, the flux derived by SNO 
using CC events with energy above 6.75 MeV should be consistent with the flux 
measured by Super-Kamiokande using ES interactions above 8.5 MeV.feil Instead, 
correctinff for the different energy threshold used in the latest Super-Kamiokande 
analysis, e3 the two fluxes differ by (0.53 ±0.17) x 10^ cm~^s~^, that is to say 3.1 a. 
Therefore the SNO data exclude the possibility that the solar neutrino deficit can 
be explained by pure Ve — > J^s oscillations. 

In the hypothesis that the observed solar neutrino deficit is due to oscillations 
of I'e's to other active neutrino species {v^ and Vr), the flux of non-electron active 
neutrinos, $^t, can be inferred. The obtained correlation between <I>^t- and $e is 
shown in Figjl^. The best flt corresponds to a non-electron active neutrino flux of: 

= (3.69 ± 1.13) X 10*^ X cTO-^s-i . (44) 

The two diagonal bands correspond to constant total flux of active neutrinos, as 
measured by SNO and Super-Kamiokande (solid line) and as predicted by the SSM 
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Fig. 10. Correlation between the non-electron active component of the solar *B flux, $^r, and the 
electron component <I>e.p|The two shaded bands represent the la bapds of the neutrino fluxes as 
obtained from SNO CCB and Super-Kamiokande ES measurements £3 The ellipses are the 68%, 
95% and 99% joint probability for the two fluxes. The two diagonal bands are the measured 
(SK+SNO) and predicted (SSM)Ej bands of constant total flux of active *B neutrinos. 



model (dashed line). The measured value, ^{vx) = (5.44 ± 0.99) x 10^ cm~^s~^, is 
in very good agreement with predictions. 

In conclusion, SNO has presented a first direct indication for a non-electron 
active component in the solar neutrino flux. Moreover the experiment measures a 

neutrino flux consistent with the SSM expectations. These are only the flrst of 
several very interesting measurements that SNO should be able to perform over the 
next few years: we now look forward to the results from the experiment's second 
phase. 

Borexino,c2l a 300 ton liquid scintillator detector with PMT readout, is being in- 
stalled at the Gran Sasso Underground Laboratory. Thanks to the high radiopurity 
and to the high light yield of the scintillator, it will be able to detect neutrino- 
electron scattering with an energy threshold as low as 250 keV. This means that 
the ''Be lines, especially that at 0.863 MeV, will be within the reach of the experi- 
ment. The monochromaticity of the beryllium line should facilitate the observation 
of VO-induced variation of the solar neutrino flux, which for Borexino would be a 
spectacular effect. Moreover Borexino will also be able to check the LOW solution, 
by studying the day/night effect for beryllium neutrinos, which is expected to be 
strong if indeed oscillations occur with the LOW parameters Jl3 Borexino is expected 
to start data-taking in 2002. 

A test of the LMA solution will be made by the KamLANDc3 experiment, a 
1 kton liquid scintillator detector, located in the Kamioka mine in Japan. The 
experiment will detect the low-energy (1 — 8 MeV) De flux emitted by the nuclear 
reactors at several Japanese power stations, at distances between 150 and 210 km 
from the detector. The estimated sensitivity of KamLAND is such that, at 90% C.L. 
and after three years of data taking, the experiment should be able to cover the entire 
domain defined by the LMA solution (down to Am^ ~ 4 x 10~^ eV'^). If oscillations 
are observed, the oscillation parameters should be determined to a precision of 20%, 
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at 99% C.L.. KamLAND is expected to start data-taking at the end of this year. 

More experiments are planned to study solar neutrinos in the future. Among 
them one should certainly mention ICARUSJ13 a liquid-argon TPC detector, which 
will be mainly sensitive tci|_£B and of which the first 600 tons half-module is cur- 
rently under test, HERON which plans to study pp and beryllium neutrinos using 
superfluid helium as a target, and HELLAZ,e3 which intends to study low-energy 
neutrinos, down to the pp neutrino energies, using a gaseous helium TPC. 

In the previous sections we have outlined the solar neutrino problem, which his- 
torically represents the first experimental hint for neutrino oscillations. We have 
described the standard solar model results and the arguments which make the model 
convincing and robust. We have discussed how experimental data are in striking 
contradiction with expectations, if no new physics is taken into account, and we 
have summarize the results obtained when the solar neutrino deficit is interpreted 
as being due to neutrino oscillations. We now want to plunge into another exit- 
ing domain, that of atmospheric neutrinos, which undoubtedly provide the most 
convincing experimental evidence for neutrino oscillations. 



6. The Atmospheric Neutrino Anomaly 



6.1. Atmospheric Neutrino Flux and Neutrino Oscillations 

Atmospheric neutrinos are produced in the decay of secondary particles created 
in the interaction of primary cosmic rays-' with the Earth's atmosphere. If the energy 
of the secondary particles is sufficiently low (< 2 GeV) that all of them decay, we 
have: 

p+Af >TT^+X 

7r± — + (45) 

Assuming that no effect can change the flavour composition of the shower before it 
is measured on Earth, Eq.(45) implies that: 

^-^^^^-2 (46) 

where AC,, (AC^ ) A/p,, (A/p^ ) are the number of niuon (electron) neutrinos and 
anti-neutrinos respectively. The exact value of TZ can in principle be affected by 
several effects, such as the primary spectrum composition, the geomagnetic cut-off, 
solar activity and, of course, the details of the model for the development of the 
hadronic shower. However, it has to be said that, although the absolute neutrino 
fluxes are rather badly known (predictions from different calculations disagree by 
~ 20 — 30%), the ratio ( ^6| ) is to ffist order insensitive to such uncertainties and is 

^ mainly protons (^^ 80%) and o-particles (~ 15%), plus a small contribution from heavier nuclei .H 
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known to ~ 5%. Neutrino oscillations would manifest themselves as a discrepancy 
between the measured and the expected value of the ratio 7?,. 

In order to quote a number which is independent of experimental parameters, 
such as energy thresholds and detector acceptances for signal and background, the 
result is generally presented in terms of the double ratio: 



7^' = 









Ti-DATA 




DATA 















(47) 



where TZdata is the measured ratio of muon-like (i.e. v^, D^) over electron-like 
(i.e. i>e, i^e) events and TZmc is the same ratio as obtained from a Monte Carlo 
simulation. Assuming a correct modelling for TZmc, Ti' <1 could either mean that 
there is a deficit of measured muon-like events, or that an excess of electron-like 
events has been observed. The quantity (|4^) alone cannot discriminate between 
these two possibilities. Moreover, a low value of TV is not by itself a proof for neu- 
trino oscillations, since one could imagine other mechanisms (i.e. neutrino decay) 
inducing a similar effect. 




Fig. 11. Sketch of the atmospheric neutrinos detection principles. 

Another, more sensitive way to detect neutrino oscillations is the study of the 
zenith angle distribution of the incoming atmospheric neutrino. Downward-going 
neutrinos, produced in the atmosphere above the detector, will travel a path of the 
order of 10 — 20 km, while upward-going neutrinos, produced at the opposite side of 
the globe, will have travelled up to 12, 000 km before detection. This means that, 
by studying the neutrino flux as a function of the zenith angle, one has access to 
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baselines spanning three orders of magnitude and can therefore hope to observe the 
modulation induced on the flux by neutrino oscillations. 

If no oscillation occurs between the production point and the detector, assum- 
ing for simplicity that no other effect alters the cosmic ray angular distribution, the 
neutrino flux per unit area is the same at any location on the Earth. If we con- 
sider the contribution to the flux from a point at zenith angle O (downward-going 
neutrinos) and we compare it to the flux at zenith angle n — Q (upward-going neu- 
trinos), based on simple symmetry considerations (see Fig.|TT|) , we can demonstrate 
that the two quantities must have the same value in the absence of new physics. 
In summary, in the no-oscillation hypothesis, the zenith angle distribution must be 
up-down symmetric, assuming no other phenomena affecting the neutrino angular 
distribution relative to the local vertical direction. Conversely, any deviation from 
up-down symmetry could be interpreted as an indication for neutrino oscillations. 
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Fig. 12. Compilatjon of the TZ' ratio (see text) measured by several atmospheric neutrino experi- 
ments (from ref.Cj). 

Historically, the first results on what is nowadays known as the atmospheric 
neutrino anomaly came from experiments originally designed to search for proton 
decay, for which the atmospheric neutrino flux constituted a background source. A 
cornnilation of results on TZ' from different experiments is shown in Fig.|l2| (from 
ref .eZI) : all measurements are systematically below one and all but two give a value 
for 7^' close to 0.7. 

The first experiment to report, in 1986, a discrepancy between the observed 
and the predicted number of atmospheric neutrinos, although that was not_im- 
mediately recognised as a possible effect of neutrino oscillations, was IMBJIj a 
water- Cherenkov detector located in the Morton mine, Cleveland, Ohio, USA. Two 
years later Kamiokande confirmed that the measured deficit of muon-like events 
was the order of 30%.l3 However, the same effect_was not observed by two other 
proton decay experiments, NUSEXj^ and Frejus,Ell, both using fine-grained iron 
calorimeters. This discrepancy led for some time to the belief that poorly under- 
stood systematic effects, mainly related to an incomplete description of neutrino 
interactions in iron and water, were inducing an intrinsic difference between the 
two experimental techniques. Although later calculations had already shown that 
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the neglected physicalnrocesses could account for at most a 10— 15% change in the 
expected event ratesj^ the issue was only definitively resolved when the Soudan2 
experiment,E3 another fine-grained iron calorimeter, located in the Soudan mine 
(2100 m.w.e.), Ely, Minnesota, USA, confirmed Kamiokandc and 1MB results. It 
is now common belief that the original discrepancy was due to fluctuations in NU- 
SEX and Frejus data. In 1994 Kamiokandc showed a distortion of the zenith angle 
distribution in the multi-GeV contained events (see below for a definition) which 
seemed to confirm the neutrino oscillation hypothesis. Kamiokandc and 1MB as well 
as Baksan,L3 an experiment in the former Soviet Union, and MACRO jl a proton 
decay experiment at the Gran Sasso Laboratory, also tried to extract information 
from the sample of particles travelling in the vertical direction, but the results were 
not consistent with each other and none of them seemed conclusive. 

The big change in the scientific community's perception of the atmospheric neu- 
trino a.nomal y cert ainly camc from the results obtained by the Super-Kamiokande 
experiment ,e3E3'c2I which in summer 1998 announced evidence for neutrino oscilla- 
tions in the atmospheric neutrino sample. However, before we discuss atmospheric 
neutrinos in Super-Kamiokande, we shall first present some of the ideas and the 
status of the atmospheric neutrino flux calculations. 

6.2. Atmospheric Neutrino Flux Calculations 

The neutrino flux on Earth (approximately 1 event /cm^/sr/sec) can in principle 
be calculated either starting from the spectrum of the primary protons and alpha 
particles which initiate the hadronic cascade in the atmosphere, or from the mea- 
sured muon flux at high altitudes. In the first case it is necessary to parametrize 
the hadronic shower development, thus requiring a relatively good understanding 
of the pion production process at the relevant energies. On the other hand, the 
muon acceptance has the disadvantage of being quite sensitive to variations of the 
geomagnetic field. Moreover the exposure time for muon flux measureroiinlsJis quite 
short, the relevant experiments being generally performed on balloons j^HS while 
primary fluxes measurements, made at high all^^bove the atmosphere, are 
usually made over relatively longer intervals .EiKj^ESE2H2a In any case a measure- 
ment of the muon flux gives a useful constraint on the the atmospheric neutrino 
flux calculations. 

At low geomagnetic latitudes the trajectories of low energy primary cosmic rays 
can be significantly deviated under the action of the geomagnetic field and, in ex- 
treme cases (energies of the order of few GeV), they would not reach the atmosphere 
at all. At intermediate energies (10 — 20 GeV), the flux of primary cosmic particles 
exhibits a strong east-west asymmetry, while at high energies the flux is essen- 
tially unaffected by the terrestrial magnetic fleld. In summary, the net effect of 
the geomagnetic field is a non-isotropy of the neutrino flux even in the absence of 
oscillations. This has been observed by Super-KamiokandeH which has measured 
an east-west asymmetry (see Fig.|l^, from refrS) consistent with expectations, thus 
confirming that the details of the geomagnetic effects on the cosmic radiation are 
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Fig. 13. East-west asymmetry measured by Supcp-KamiokandeM as a function of the lepton 
momentum (from ref.Q), compared to expectations.Lj and are the number of events from 
the east and the west respectively. 



well understood. In Super-Kamiokande (geomagnetical latitude 25.8'^ N, close to 
the geomagnetic equator) the cutoff momentum for protons arriving horizontally 
from the east is ^ 50 GeV. 

A modulation of the low energy component of the primary cosmic ray flux is 
expected as a function of the solar activity (cycle of ~ II years), which induces 
fluctuations in the geomagnetic field (see for example ref.EZl): the greater the solar 
activity, the lower the primary cosmic ray flux on the Earth. The correlation be- 
tween the relative sun-spot number and the cosmic rav flux can be studied using 
the data from the Mount Washington neutron monitor .E2 The influence of the solar 
wind is only sizeable (~ 15 — 20%) for neutrino energies of 1 GeV or below and 
only at high geomagnetic latitudes, where the geomagnetic effect is minimal, while 
it becomes negligible at ~ 2 GeVo 

The experimental results for atmospheric neutrinos are usually presented in 
terms of the neutrino flux calculations performed by Honda and collaborators 
(HKKM)EJ and by the Bartol groupO A comparison between the two models and a 
compilation of some of the most recent data on primary protons is shown in Fig.p^ 
(from refJIj). 

Both HKKM and the Bartol models are based on one-dimensional computations, 
in the sense that the momenta of all secondary particles, hadrons, charged leptons 
as well as neutrinos, are assumed to be collinear with that of primary cosmic ray 
(1-D calculation). It has only very recently been realized that three-dimensional 
(3-D) calculations can be important for a correct description of the neutrino spatial 
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distribution, especially in the sub-GeV energy rangeOo'ES The new calculations, 
originally motivated by the need of having a better description of the hadronic pro- 
cesses involved in the shower development, are based on the the FLUKA simulation 
package. L3 The effect of taking the three-dimensional development of the shower 
at low energies (< 1 GeV) can be clearly seen in Fig.|l^, where the results from 
FLUKA-based 1-D and 3-D simulatinns are shown for the zenith angle distribution 
at the Kamioka site (figure from refJid). 

The Bartol and FLUKA groups have started a comparison of the two models, 
using the same input primary spectrumE3 The results of this study show that 
the two normalisations disagree by approximately 20% (Bartol is lower) and the 
difference between the two models appears to be energy-dependent. According to 
the author of ref.Ol, a preliminary comparison with HKKM shows that, if they 
had used the same primary spectrum as an input, also their predictions on the 
normalisation would have been lower than FLUKA's. 

In summary, even if the e/fj, ratio is rather well known, the uncertainties on 
the absolute neutrino fluxes are still quite large. Among the dominant sources 
of systematic errors are the knowledge of the primary fluxes, still limited despite 
new data has been recently made available, and the uncertainties on the secondary 
particle production model. Other effects, such as the atmosphere modelling and 
the detector altitude, play only a marginal role. In order to improve the theoretical 
calculations it is very important that further constraints on the models are obtained 
from new experimental data. 

Several experiments relevant for the atmospheric neutrino flux calculations should 
be taking data over the next few years. HARP,lZI an experiment running at the 
CERN PS, aims to study hadroproduction in a range of energy relevant for future 
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Fig. 15. Comparison between 1-D and 3-D predictions for the neutrino fluxes at the Kamioka site, 
as a function of the zenith angle (from ref.M). 



projects such as neutrino factories and for the determination of the atmospheric 
neutrino flux. More recently a proposal for another experiment has been submitted 
to CERnS if approved, that experiment would make use of the NA450 appara- 
tus to study hadroproduction in a range of energies (120 GeV) relevant for both 
atmospheric neutrinos and long-baseline neutrino oscillation projects. 

Among the new proposals to study cosmic ray muons for the determination of 
the neutrino fluxes, an interesting project, known as adlerB is now in its design 
stage. The main idea is to fly a small detector on board an aircraft to measure 
the atmospheric muon flux at different geomagnetic latitudes. This possibility, 
together with the relatively longer exposure times, would be an obvious advantage 
with respect to balloon-borne experiments. 

6.3. Atmospheric Neutrinos in Super- Kamiokande 

We have already given a brief description of the Super-Kamiokande experiment 
and will not repeat it here, however it is worth summarising some information 
relevant for the following discussion. Atmospheric neutrinos are detected in Super- 
Kamiokande by measuring the Cherenkov rings generated by the primary particles 
produced in the neutrino CC interactions with the water nuclei. Thanks to the high 
PMT coverage, the experiment is characterised by an extremely good light yield (~ 8 
photo-electron per MeV) and can detect events of energies as low as ~ 5 MeV. The 
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large detector mass and the possibility of clearly defining a large inner volume allow 
to collect a high statistics sample of fully contained events (FC) of relatively high 
energies (up to ~ 5 GeV), the FC events being defined as those having both the 
neutrino vertex and the resulting particle tracks entirely within the fiducial volume. 
The contamination from downward-going cosmic muons is drastically reduced by the 
containment requirement on the primary vertex coordinates. Fully contained events 
can be further subdivided into two subsets, the so-called sub- GeV and multi-GeV 
events, with energies below and above 1.33 GeV respectively. In Super-Kamiokande 
jargon FC events include only single-ring events, while multi-ring ones (MRING) 
are treated as a separate category. Another sub-sample, defined as the partially 
contained events (PC), is represented by those charged-current interactions where 
the vertex is still within the fiducial volume, but at least a primary particle, typically 
the muon, exits the detector without releasing all of its energy: for those events the 
energy resolution is therefore worse than for FC interactions. Finally, upward- 
going muons (UPMU), produced by neutrinos coming from below and interacting 
in the rock, can also be used to independently check the neutrino oscillation result. 
In the literature they are sometimes further subdivided into stopping muons and 
through-going muons, according to whether or not they stop in the detector. The 
different samples defined aboiie explore different ranges of the neutrino energy: this 
is shown in Figjl^ (from ref.E3), where the event rates for sub-GeV and multi-GeV 
FC events, as well as for through-going and stopping muons is shown. For PC 
events the neutrino energy is in the multi-GeV range. 




Fig. 16. Event rates as a function of the neutrino enp«y for some of the different categories of 
neutrino interactions defined in the text (fig. from refE2l)- 



Particle identification in Super-Kamiokande is performed using likelihood func- 
tions to parametrize the sharpness of the Cherenkov rings, which are more diffused 
for electrons than for muons. The algorithms, which h ave been tested on cosmic 
muons, decay electron samples and using test-beam data,Eil are able to discriminate 
the two fiavours with very high purity (of the order of 98% for single track events). 

The most recent value for TZ' reported by the Super-Kamiokande collaboration, 
based on 1289 days of data, is O.GSS^o oiy ± 0.050 for the sub-GeV sample and 
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0.675t[J;J^^2 ± 0.080 for the multi-GeV sample (both FC and PC).!! Thus the mea- 
sured value of TZ' is different from unity by about 7 a. 
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Fig. 17. Distribution of the cosine of the zenith angle obtained by SuperJ4amiokade from 1289 live 
days data, for electron-like and muon-like contained events (from refo). The solid line is the 
distribution expected in absence of oscillations, while the hatched histogram is for i/^ — > i/^ 
oscillations in the two-flavour mixing scheme, with maximal mixing and sin 26 = 2.5 X 10~3 eV^. 
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The distributions of the cosine of the zenith angle for the 
multi-GeV samples are shown in Fig.^ (1289 live days, from ref 
sets the electron-like events are in good agreement with predictions in absence of 
oscillations, while the muon spectrum, which strikingly disagrees with no-oscillation 
expectations, is well described by an oscillation induced modulation. It has to be 
noticed that, while the sub-GeV muon sample is depleted over the entire spectrum, 
the deficit being larger for smaller values of cosO, the spectrum of neutrinos com- 
ing from above (cos > 0) is very little changed in the case of multi-GeV muons. 
This is consistent with oscillations occurring with a value of Am^ which is better 
matched by low-energy neutrinos (sub-GeV), while for high energy neutrinos the ef- 
fect of oscillations becomes sizeable only for larger baselines (neutrinos coming from 
below) . However it has to be stressed that the angular resolution degrades at lower 
energies, where the lepton momentum direction becomes a worse approximation for 
the incoming neutrino direction. 



A useful dimensionless quantity which can be derived from the zenith angle 
distribution is the so-called up-down asymmetry A, defined as the difference between 
the number of upward-going (U, cosO < —0.2) and the number of downward-going 



36 An Experimentalist's View of 



(D, COS > 0.2) neutrinos, divided by the sum of the two numbers: 

We have seen that, neglecting the geomagnetic effect corrections, which can be taken 
into account in the simulation, there is no reason to have any asymmetry in the 
neutrino flux. Any deviation of A from zero is a signature for new physics and can 
be interpreted in terms of neutrino oscillations. The dependence of this quantity 
on the measured particle momentum is shown in Fig.^ (left) for electron-like and 
muon-like events (both FC and PC). The uncertainty assumed on the absolute 
normalisation is 25%, with 20% contribution from the neutrino flux calculation and 
15% from the neutrino interaction cross section. Also here the experimental data are 
consistent with neutrino oscillations in the muon sample, while electron neutrinos 
appear to be unaffected. In agreement with the oscillation hypothesis, we see that 
for muon-like events the asymmetry is essentially zero at low energies, where the 
oscillation length is small and therefore both upward-going and downward-going 
muons are equally depleted, while for increasing values of the muon momentum, 
when mostly neutrinos from below are affected by the oscillation, the asymmetry 
becomes more and more negative. 
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Fig. 18. Super-Kamiokande (1144 days) up-down asymmetry for electron-like and muon-like 
events, as a function of the particle momentum(left) and ratio between observed and expected 
number of electron-like and muon-like events as a function of L/E (right). The dashed lines corre- 
spond to U/j, — > Ut oscillations with maximal mixing and Am^, = 3.2 X 10~^ eV^, while the hatched 
boxes indicate no-oscillation expectations (figures from refO). 



Another way of looking at the same effect is to consider the ratio between ob- 
served and expected number of electron-like and muon-like events as a function of 
the ratio L/E: this is shown in Fig.|l^ (right) for 1144 live day of Super-Kamiokande 
data. Also here electron-like events show no evidence for oscillations, while muon- 
like events show a dependence on L/E which, although the oscillatory pattern is 
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washed out by angular and energy resolution, is consistent with the oscillations 
hypothesis. 

Interpreting the muon-like event deficit as the result of — > Vr oscillations in 
the two-flavour mixing scheme, Super-Kamiokande extract an allowed domain for 
the oscillation parameters. Events are binned in a multi-dimensional space defined 
by particle type, energy and zenith angle, plus a set of parameters to account 
for systematic uncertainties. The result of the fit based on 1289 days of data. 



eft. The best fit 
For comparison, 



using FC, PC, UPMU and MRING events, is given in Fig.|l 
corresponds to maximal mixing and Am^ = 2.5 x 10~^ eV^.Ei 
an earlier plot*"' of Super-Kamiokande_allowed reeiDn (90% C.L.) is^shown together 
with those obtained by Kamiokande,E3, Soudan2,E3 and MACRO .E3 In particular it 
is interesting to consider how the preferred solution for Am^ indicated by Super- 
Kamiokande is lower than that previously presented by Kamiokande, with only a 
very small overlap between the 90% C.L. contours obtained by the two experiments. 
The reason for this effect is not fully understood, although one possible explanation 
could be the different fitting procedure, which for Kamiokande is based on the 
flavour ratio TZ' only, while in Super-Kamiokade a shape analysis of the energy and 
the zenith angle spectra is performed. 

The Soudan2@ and the MACRoEeI experiments, although with less sensitivity, 
have confirmed an atmospheric neutrino anomaly consistent with that observed by 
Super-Kamiokande . 
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Fig. 19. (Left) Allowed region at different C.L. a«-«btained by Super-Kamiokande (1289 days) 
in the two-flevour i^^ — > Vt hypothesis (fp«m refE3).i— /Right) An e»#jier contour from SuBer- 
KamiokandeEZI is compared to Kamiokande,Cj, Soudan2E3 and MACROEj results (fig. from ref.Lj). 



The Super-Kamiokande interpretation of the atmospheric neutrino anomaly in 
terms of muon disappearance rather than electron appearance is corroborated by 

*the arguments outlined in the following do not change if the latest Super-Kamiokande results are 
considered. 
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the results obtained by CHOOZEI and Palo Verdej^a two neutrino oscillation experi- 
ments at nuclear reactors, which have studied 9^ disappearance in the parameter re- 
gion relevant for atmospheric neutrinos. In particular CHOOZ limits (at 90% C.L., 
Am^ < 7 X 10^^ for maximum mixing and sin^ 29 < 0.10 for large Am^) seem to 
exclude i'^ <-> Ve mixing in the region of interest, at least as a dominant modeE3 
However the ambiguity between oscillations to tau or sterile neutrinos is not resolved 
by the analysis of Super-Kamiokande data described previously. In order to discrim- 
inate between the — > Vr and the i'^ ~f Ug hypotheses, it is necessary to consider 
other physics effects able to distinguish one mode from the other. First of all, while 
the v-T has standard neutral current interactions with matter, does not couple to 
the boson, with a consequent depletion of NC events and an up-down asymmetry 
in the NC sample. Moreover, in case of mixing with a sterile neutrino, matter effects 
- see Eqs.(p6|), ( p7| ) and ( p8|) - would induce a sizeable change in the.oscillation prob- 
ability for neutrino energies such that E/\/S.'m?\ > lO'^GeV/eV^.Eil Based on those 
considerations, Super-Kamiokande has published an analysis of 1100 live days data, 
showing that the u^^ Vt hypothesis fits loarenergy CC interactions, but does not 
fit the NC sample or high energy CC data.tS Although they cannot exclude more 
complicated schemes in which the mixing is between i'^ and admixture of v-^ with 
a small component of pure i/^ is disfavoured by their data at 99% C.L.. 

Another way of testing the v^^ — > Vs hypothesis, which is again motivated by the 
not having standard electro-weak interactions, is to look for possible suppressions 
of tt" production. This technique is currently limited by the systematic error on the 
TT*^ production cross-section^ but the situation should improve in the near future, 
since the K2K experiment j2j which we shall discuss shortly, is expected to reduce 
the uncertainty from the present value of ~ 20% down to ^ 5%. 

Super-Kamiokande has recently presented preliminary results on Vt appearancc.E 
The experiment does not have the resolution to unambiguously identify an event as 
being a Vt CC interaction, however they have developed a likelihood-based analysis 
to extract a sample of tau-like events. The observed number of selected events is 
419, for an expected number of 31 j^t's plus 387 background events; moreover the 
measured zenith angle distribution is better fitted by 's plus background than by 
background only. 

In summary, we have seen that Super-Kamiokande results, which have been 
confirmed with less sensitivity by other atmospheric neutrino experiments, represent 
a strong hint for neutrino oscillations with near-to-maximal mixing (sin^ 2Q >_p.88) 
and Am^ in the range 2 — 6 x 10"'^ eV^. New experiments, like ICARUSe3 and 
MONOLITHo will study or plan to study atmospheric neutrinos in the future. 
However the uncertainties on the atmospheric neutrino flux predictions are such 
that a firmer experimental proof of neutrino oscillations is required, together with a 
better determination of the oscillation parameters if the effect is confirmed. For this 
reason it is very important to check the atmospheric neutrino anomaly in controlled 
experimental conditions, by using artificially produced neutrinos (i.e. accelerator 
beams) . 
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Typical neutrino energies at accelerators are in the range ^ 1—30 GeV, therefore, 
in order to match the Am^ domain suggested by atmospheric neutrino results, it 
is necessary to increase the detector-to-source separation to distances of the order 
of hundreds of kilometres. Experiments of this sort are commonly referred to as 
long-baseline experiments and we shall discuss them in the next section. 

6.4. Long-baseline Accelerator Experiments 

The first accelerator-based long-baseline neutrino oscillation experiment ever 
built is K2K (KEK-to-Kamioka) jHj which makes use of an almost pure beam 
(98.2%) directed from KEK, in Japan, to the Super-Kamiokande detector, 250 km 
away. K2K, which has been taking data since June 1999, is a disappearance experi- 
ment, aiming to measure a deficit of muon neutrinos in the far detector at Kamioka 
compared to the initial beam intensity measured in the near detector at KEK. The 
neutrino beam, produced in the interactions of 12 GeV primary protons from the 
KEK-PS on an aluminium target, has an average energy of ~ 1.3 GeV (hence the 
maximum sensitivity to oscillations is for Am^ ^ 5 x 10^"^ eV^). Pion and muon 
detectors, placed along the beamline, are used to measure the secondary particles 
after the horn focusing system and to monitor the beam centering and its inten- 
sity. The design proton intensity is 6 x 10^^ protons per pulse (1.1 spill/2.2 s 
period). The results recently presented by the K2K collaboration are based on 
2.29 X 10^^ pot (protons-on-target), collected between June 1999 and June 2000. 
The alignment between the neutrino target and the far detector is performed using 
a GPS survey, giving a precision better than 0.01 mrad, while the precision of the 
civil construction is better than 0.1 mrad. GPS timing information is also used to 
correlate events in Super-Kamiokande with the KEK-PS spill. 

In order to have a reliable measurement of the beam intensity at production 
and to reduce the systematic uncertainties related to the experimental technique 
adopted, the near detector has been built using the same principles and technol- 
ogy as for Super-Kamiokande. Located 300 m from the production point, it is a 
1 kton water-Cherenkov detector with photo- multiplier readout, the PMTs being 
arranged so to provide the same coverage as for Super-Kamiokande (40%). The 
water-Cherenkov detector is complemented by a scintillating fibre tracker (with a 
6 ton water target), a lead-glass electromagnetic calorimeter and a muon range 
detector. The combined information from the different sub-detectors allows mea- 
surement of the beam profile, its energy distribution and intensity, as well as the 
determination the contamination (~ 1.3%) and the study of neutrino interactions 
in the 1 GeV region. 

The beam has been measured to be centred within systematic uncertainties 
(0.7 mrad) and the centre has been shown to be stable within ±1 mrad. The 
muon energy and angular distributions are also monitored on a daily basis and so 
is the beam intensity: they all appear to be stable over the period of data taking 
considered. Given this stability, it is possible to predict the number of expected 
events at the far site in absence of oscillations by extrapolating the normalisation 
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measured at the near site to the far detector location. The beam simulation used 
for the extrapolation is validated by a comparison between the predicted secondary 
pion momentum and angular distributions and those measured at the pion monitor 
placed along the beamline. The selection criteria used to identify muon-like and 
electron-like interactions at the near detector are essentially the same as those used 
in Super-Kamiokande for the FC atmospheric neutrino analysis, with the addition 
of a cut on the total number of photoelectrons to discard spills with multiple events. 

Table 5. Number of expected and observed events in the far detector of the K2K experiment. 



Sample 


Expected Events (No osc.) 


Observed Events 


1-ring fi-lihe 


20.9 


14 


1-ring e-like 


2.0 


1 


MRING 


14.9 


13 


Total 


37.8 ± 0.2(stat.)ti%isys.) 


28 




Fig. 20. K2K expected sensitivity for 10^'' pot, compared with the allowed regions irf4tie oscillation 
parameter space obtained by Super-Kamiokande, Soudan2 and MACRO (from ref.ESI). 

The systematic uncertainty on the near detector flux measurement is 5%, with a 
statistical error of < 1%. For the far detector the systematic error is 3%, while the 
uncertainty related to the near-to-far extrapolation is estimated to be 1^%. Other 
uncertainties, such as those on the beam spectrum and the neutrino interaction cross 
sections, are negligible to first order, due to cancellations in the near-to-far compari- 
son. The expected number of events in absence of oscillations is 37.8±0.2{stat.)~^f^, 
while only 28 have been observed in real data. In Tab.|| the predicted and observed 
number of events are compared for the different categories (one-ring muon-like, 
one-ring electron like and multi-ring) : while essentially no conclusion can be drawn 
about electron-like events, due to the small statistics, the niuon neutrino sample 
appears to be significantly depleted, a fact which might be interpreted as evidence 
for neutrino oscillations. However, the K2K collaboration has not yet given any in- 



An Experimentalist's View of . . . 41 



terpretation of their result and is waiting to have more statistics (10^" pot by year 
2005) and a better understanding of systematic uncertainties before drawing any 
strong conclusion from the observed anomaly. The expected final sensitivity for K2K 
is shown in Fig.^, together with the allowed solutions obtained by Kamiokande, 
Super-Kamiokande, MACRO and Soudan2. From that plot it becomes clear that, 
unless the actual Am^ is at the high end of Super-Kamiokande allowed region, K2K 
will not provide a definitive proof that the atmospheric neutrino anomaly is due to 
neutrino oscillations. 

Two other long-baseline projects are expected to start operation over the next 
few years: one in the USA, the MINOSEj experiment at Soudan, using the NuMEJ 
beam fromFermilab, and the other in Europe, using the CERN-to-Gran Sasso neit 
trino beamO (NCGS) in c onju nction with the two proposed experiments OPERAEj 
and ICARUS/ICANOE.S'H 

The MINOS experiment aims to investigate the atmospheric neutrino anomaly 
by measuring the neutrino beam at two different locations, one close to production 
(near detector, at Fermilab, ~ 320 m from the production point), where no effect is 
expected from neutrino mixing, and the other far away (far detector, in the Soudan 
mine, Ely, Minnesota, about 730 km from Fermilab), where neutrino oscillations 
would modify the beam flavour composition and induce distortions in the measured 
energy spectra. 




Fig. 21. Predicted event rates at the MINOS far site for the three possible beam configurations, 
together with the spectrum expected for "perfect focusing" , where all produced neutrinos reach 
the detector. 

The NuMI neutrino beam is obtained from the decay of secondary particles 
produced in the interactions of 120 GeV protons, delivered by the new Fermilab 
Main Injector, with a carbon target. The expected proton flux is 4 x 10^'^ protons 
per pulse, with repetition rate 1.9 s and spill duration time 10 /iS. The secondary 
mesons are focused by two parabolic horns and allowed to decay in a 675 m long 
decay tunnel. The relative distance between the two magnetic elements can be 
changed by moving the second horn with respect to the other: by adjusting the 
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horn positions and currents it is then possible to select different regions of the phase 
space, thus providing the flexibility to optimize the beam energy during the course 
of the experiment. Three beam settings arc possible (see Fig.^) , corresponding to 
three different energy regimes: low (< i?^ >~ 3 GeV), medium (< E„ >~ 7 GeV) 
and high (< E^, >~ 15 GeV) energy. In order to have optimal sensitivity in the 
low part of the atmospheric neutrinos allowed region, which seems to be preferred 
by Super-Kamiokande data, MINOS has decided to start taking data with the low- 
energy beam configuration. Another focusing system is under consideration and 
might be included in the final design of the NuMI beam, the so-called hadron hose. 
The basic idea is to have a wire carrying a 1 kA current down the centre of the 
decay tunnel, to provide additional focusing for the charged secondary mesons, with 
the double advantage of increasing the beam intensity and reducing the systematic 
uncertainties on the beam extrapolation from the near to the far site. 

In order to minimize the impact of systematic uncertainties, MINOS far and 
near detectors (5.4 kton and 980 ton mass respectively) have been designed to be as 
similar as possible. They are both magnetised iron-scintillator sampling calorime- 
ters with photo- multiplier readout. The average toroidal magnetic field in the iron 
is about 1.3 T. The longitudinal and transverse granularity is the same for the 
two detectors: 2.54 cm thick octagonal iron plates are alternated with scintillator 
planes, made of plastic slabs with 1x4 cm^ transverse section and lengths varying 
between 4 m and 8 m for the far detector planes. The slabs orientation on successive 
planes differs by OO*', so to provide two coordinate measurements on the transverse 
plane. A green wave-length shifting fibre, 1.2 mm diameter, is glued into a grove 
at the centre of the wide side of each slab: the light produced in the scintillator is 
trapped into the fibre, shifted in frequency and then routed to multi-pixel PMTs' 
through a clear fibre, optically connected to the wave-length shifting one. The read- 
out is two-sided at the far detector and one-sided at the near detector. While the 
far detector planes are all fully equipped with scintillator strips, the near detector 
will be only partially instrumented. In the forward part (the first 120 plates, out 
of a total of 180), which logically defines the target and calorimeter sections, only 
every fifth plane is fully instrumented, while 4/5 of the plates are instrumented in 
a limited region around the beam fiducial area. In the remaining section, the muon 
spectrometer, every fifth plane is fully instrumented. 

A calibration detector, much smaller in size but with the same longitudinal 
segmentation as the MINOS far and near detectors, is being tested at the CERN 
PS,E£il to characterize its response to muons, hadrons and electrons and to study 
the performances of the different hardware components used by MINOS. 

The main physics goals of MINOS, which is expected to start taking data in 
2004, are to test the neutrino oscillations hypothesis, to verify that the dominant 
mode is — > as suggested by Super-Kamiokande data, and to measure the 

' 16-pixel PMTs are used at the far detector, where the signal is optically multiplexed, while at 
the near detector 64-pixel PMTs have been chosen and multiplexing is applied only in the muon 
spectrometer section. 
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oscillation parameters with a precision of about 10%. Moreover the i/^ 
i/^ — > Uf, channels will also be investigated. 



and 



The criterion envisaged by MINOS to verify the oscillations hypothesis is the 



so-called T-test, based on the parameter 



T = 



(49) 



where A/i^ and A/q^ are the number of muon-like and muon-less events respectively. 
If J'^'s oscillate to another active neutrino species, say i^r's, the muon-like sample 
will be depleted {disappearance of v^) and the measured A/i^ in the far detector 
will be less than that expected in the absence of oscillations; on the other hand, 
VtGG events will look very much like neutral current interactions and, for neutrino 
energies above the r production threshold, will populate the muon-less sample, 
so that the total number of events will remain unchanged. Therefore, — > v-r 
oscillations would manifest themselves as a low Tfar/Tnear ratio with respect to 
expectations in absence of oscillations. If the I'^'s oscillated to sterile neutrinos, the 
muon-like and the muon-less samples would be depleted in the same manner, since 
sterile neutrinos do not couple with any of the weak interactions mediating bosons. 
Therefore for pure i'^ — > Vs oscillations the ratio Tpab/Tnear would be the same 
as for the no-oscillations case. However, the NC sample would also be affected by 
oscillations to sterile neutrinos, which would induce a reduction in the number of 
neutral current interactions and a distortion of the NC energy spectrum: if none of 



those effects is measurable, a limit can be set on the Vu 



oscillation amplitude. 



CC energy distributions — Ph2le, 1 kt.yr,, sin^(2i?) = 0,9 

Am'=0.002 eV' Am'=0.0035 eV' Am^=0.005 eV' 
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Fig. 22. (Top) Energy distributions at the MINOS far detector, in absence of oscillations (solid 
line) and as would be measured for different values of Am^ and sin^ 2d = 0.9. (Bottom) The 
ratio between the oscillated and the unoscillated spectra shows the characteristic dip shape, at an 
energy determined by the value of Am'^ and whose depth depends on the oscillation amplitude. 

If the oscillation hypothesis is confirmed by the T-test, assuming the dominant 
oscillation mode is — *■ the oscillation parameters can be extracted by com- 
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paring the v^GC energy spectrum measured in the far detector to expectations in 
absence osciUations as obtained from the measured spectrum in the near detector. 
When considering the ratio between the two distributions (see Fig.E2E2j), a very 
characteristic dip appears, at an energy corresponding to the oscillation Am^ and 
whose depth is determined by the oscillation amplitude sin^ 29. 




Fig. 23. (Left) MINOS sensitivity at 90% C.L. for Vr as obtained from the T-test; results 

are shown for lOkton X year and all three possible beam configurations. (Right) 90% C.L. allowed 
regions from fits to reconstructed CC energy distribution, for sin^ 29 = 0.9 and different values of 
Amp; results are shown for the low-energy beam only. 

If no oscillation signal is found, the 90% C.L. limits which MINOS would ex- 
tract from the.^ T-test analysis, based on 10 kton x year of data, are shown in 
Fig.^ (Ieft),ll23 together with the allowed regions obtained by Kamiokande and 
Super-Kamiokande. Using the low-energy beam, then, MINOS will cover the entire 
parameter space suggested by the atmospheric neutrino results. On the other hand, 
if an oscillation signal is observed, a measurement of the oscillation parameters can 
be performed. MINOS sensitivity to the oscillation parameters, as obtained from 
the far-to-near comparison of the i^^CC energy spectrum for 10 kton x year, is shown 
in Fig.|2^ (right) jl23 for sin^ 29 = 0.9 and Am^ varying between 0.2 x 10"'^ eV^ and 
0.5 x 10""^ eV^: as expected, the higher the value of Am^, the better the precision 
to which the parameters can be determined. 

MINOS also plan to study atmospheric neutrinos and preliminary results, based 
on a 18 kton x year sample of simulated events, show that the experiment has the 
sensitivity to check the atmospheric neutrino anomaly in the range of parameters 
suggested by Super-Kamiokande, with completely different systematics.E£il 

Intense efforts have been recently made within the MINOS collaboration to 
optimize the experiment capabilities for detection of a j/^ — > i^e signal. The main 
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difference with respect to the 



search is that in case of osciUations to Vf, an 



appearance analysis has to be performed, looking for an excess of events in the far 
detector, above the expected background from prompt ly^CC interactions and NC 
events with a topology similar to that of i>eCC. The prompt contamination in the 
beam is quite small, of the order of 0.6 — 2%, depending on beam configuration. 
The systematic uncertainties on background predictions in the far detector should 
be drastically reduced by the measurements of the same background events in the 
near detector. Moreover results from hadron-production experiment sho uld help to 



reduce the systematic error on beam predictions. Preliminary resultsE23 show that, 



if the total error can be kept within 10%, MINOS should be able to improve the 
CHOOzS limit on {Uesl^ by about a factor two.™ 

The CERN-to-Gran Sasso program is based on a completely different strategy 
than that of NuMI-to-Soudan. Instead of measuring the disappearance of muon 
neutrinos at the far location and extracting the oscillation parameters from the 
characteristic pattern induced by oscillations in the neutrino energy distribution, 
they have decided to go for direct detection (appearance) of tau neutrinos in an 
almost pure t'^ beam. The two approaches should be regarded as complementary 
and equally valuable, especially in a field where redundancy is a requirement rather 
than a choice. 
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Fig. 24. CNGS i/^ flux at the Gran Sasso site. 



CNGS is a wide-band neutrino beam of relatively high energy (< E{Uf^CC) >~ 
20 GeV, see Fig.p^ , which has been optimised to maximize the number of VrCG 
events for appearance experiments at Gran Sasso. It is expected to be commissioned 
by May 2005. The beam design benefits from the experience gained at CERN with 
the West Area Neutrino Bearn_(WANF), in most recent times used by the short- 
baseline experiments NOMAD0 and CHORUS .Ell Primary protons of 400 GeV 



are extracted from the CERN-SPS and let interact with a graphite target, where 
secondary particles, mostly pions and kaons, are produced. High energy positively 



"in the three-generation model, for |Am|g| ^ lAmJjIi ™3 > "^2 and |(7fi3p = \U-ra 
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charged mesons are focused by means of two magnetic elements, the horn and the 
reflector, and then decay in a 900 m long tunnel. The expected proton intensity 
is of 4.5 X 10^^ pot /year. The t?^, i/f, and De contamination is very low (in terms 
of CC interactions: ^ 2%, ^ 0.8% and 0.05% respectively), while the prompt i^r 
contamination (0(10^^)) is essentially negligible for the purposes of i^t appearance 
experiments. Since the baseline is fixed (about 730 km between the production 
point at CERN and the Gran Sasso Laboratory), the low value of the Am^ sug- 
gested by atmospheric neutrino experiments would require operation at the lowest 
energy possible; on the other hand, neutrinos must be energetic enough to produce 
charged tau leptons in charged current interactions. A compromise between these 
two competing requirements is achieved by adjusting the settings of the beam fo- 
cusing system. The expected nu mbe r of j/t-CC events at Gran Sasso per kton per 
year are given in TabJ^ (from ref.EHj). 

Table 6. Expected number of Ur charged current interactions for an isoscalar target 
for sin^ 29 = 1 and different values of Am^, assuming 4.5 X 10^^ pot/year (from ref. 
efficiencies have not been taken into account. 



Am2(eV2) 


N^-^CC 


/kton/year 




E„ = l- 30 GeV 


= 1- 100 GeV 


1 X 10--^ 


2.34 


2.48 


3 X 10~3 


20.7 


21.4 


5 X 10-3 


55.9 


57.7 


1 X 10^2 


195 


202 



OPERAB is an emulsion-based experiment, which aims to detect z^rCC events 
by identifying the r produced through its decay kink. Emulsions are used as high 
precision trackers rather than as the active target. The OPERA design, which 
largely benefits from the experience acquired with the CHORUsIl^ and DONUlH 
experiments, is an evolution of the Emulsion Cloud Chamber (ECC) technique, 
where the high-precision tracking capabilities of emulsions are integrated with a 
passive material to provide a large target mass. In the OPERA baseline design, an 
ECC cell consists of a 1 mm thick lead plate followed by a thin film, made of two 
50 /im thick emulsion layers separated by a 200 fim plastic base. The basic detector 
unit is obtained by stacking 56 cells to form a compact brick (10.2 x 12.7 x cm^ 
transverse section, 10 Xq length, 8.3 kg weight), which can in turn be assembled 
into walls. An additional emulsion film is placed in front of the first lead plate, 
to improve the track matching with the upstream wall, while another film, the so- 
called Special Sheet, separated from the most downstream cell by a 2 mm plastic 
plate, is used to have good r detection efficiency also for events occurring in the 
last lead plate. After each wall, an electronic tracking detector is used to select 
the brick where the neutrino interaction occurs. In the baseline design the tracking 
detector in the target section follows very closely that of MINOS, with scintillator 
strips coupled to wave-length shifting fibres and multi- anode PMT readout, where 
the selected PMTs are the same as those used in MINOS near detector. A wall 
of bricks and two planes of electronic detectors constitute one OPERA module, 



ti7- jran oasso, 
122). Detector 
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whose structure is then repeated 24 times to form a supermodule, also including 
a downstream muon spectrometer. Each muon spectrometer consists of a dipolar 
magnet (1.55 T in the tracking region), made of two vertical walls of iron layers, 
interleaved with RPC detectors, and of drift tubes, placed in front and behind the 
magnet as well as between the two walls. Three supermodules form the full OPERA 
detector, for a total target mass of about 1.8 kton. 

The data from the electronic detectors, which will be analysed quasi-online^ are 
used to select bricks where neutrino interactions occur and which have to be removed 
by an automated system. Additional bricks might need to be removed for a complete 
reconstruction of the interesting events. Events selected as t candidates will then 
be sent to dedicated scanning stations, where further studies can be performed to 
achieve the required background suppression. The expected rate of I'^CC events is 
of about 30 events per day. 

The T decay channels considered by OPERA are listed in Tab.^, together with 
the corresponding Vt detection efficiencies and the expected number of background 
events" per year. New studies are underway to improve the experiment sensitivity 
by adding also the t ^ pv channel (23.5% B.R.) to the decay modes considered. 
The expected number of r events per year is given in Tab.^ as a function of possible 
values of Am^. 

Table 7. OPERA efficiency for the r decay channels considered and the corresponding expected 
number of background events per year. The efficiencies, which also include the branching ratio for 
each decay channel, have been obtained as weighted sums of the efficiencies for DfS and non-scaling 
processes. 



Decay mode 


e{UrCC)i%) 


NbkgdI year 


T ^ e 


3.7 


0.04 


T ^ 


2.7 


0.03 


T ^ h 


2.3 


0.05 



Table 8. Expected numbes of r events observed in OPERA per year as a function of possible 
values of Hm? (from ref.E^)- 



Am2(eV2) TV^/year 



f.5 X IQ- 


-3 


0.82 


2.5 X 10~ 




2.82 


3.2 X IQ- 


-3 


3.66 



If no signal is observed, the average upper limit at 90% C.L. which would be 
obtained by OPERA is shown in Fig.|2^ (left) for 2 and 5 years of exposure respec- 
tively. On the other hand, if events are actually observed, a measurement of 
Am^ can be performed. Assuming maximal mixing and Am? — 3.2 x 10"'^ eV^, the 
90% C.L. allowed region for the oscillation parameters as determined by OPERA 
after 5 years of data taking is shown in Fig.p^ (right). 



"mostly from charm production, large angle muon scattering and hadron re-interactions. 
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Fig. 25. (Left) Average upper limit at 90% C.L. which would be obtained by OPERA if no 
oscillation signal is detected, after 2 and 5 years of data-taking. (Right) 90% C.L. allowed region 
for the oscillation parameters as determined by OPERA after 5 years of running. 



OPERA also plan to perform a —^ i>e oscillation analysis, by looking for 
appearance above the ~ 1% intrinsic Ve contamination, and to search for neutrino 
oscillations by studying the ratio between NC and CC interactions.Ej However the 
study of the systematic uncertainties for those analyses is likely to be somehow 
limited by the lack of a near detector. 

jran Sasso laboratories to 
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technology. ICARUS is a 



The other experiment which is being planned at 
detect neutrinos from NCGS is based on the ICARUS 
liquid- Argon TPC which allows event reconstruction in three dimensions, combining 
high spatial resolution with excellent particle identification capabilities and precise 
homogeneous calorimetry. ICARUS has been conceived as a multi-purpose experi- 
ment with a very wide physics program, which goes from proton decay searches to 
the study of solar and atmospheric neutrinos as well as of supernova neutrinos and, 
of course, of artificially produced neutrinos from CERN. 

The final phase of ICARUS foresees a detector with a sensitive mass of at least 
5 kton, but a step-wise strategy has been adopted to demonstrate the technical 
feasibility of such a challenging project. The first large-scale prototype had a volume 
of l,000m'^ (15 ton) and was built as a test-bench for the cryogenic system, the 
internal detector mechanics and the liquid- Argon purification. The 15 ton prototype 
was successfully tested between 1997 and 1999. The next step towards the full- 
scale detector is the T600 modulejlij which is now being commissioned at Pavia 
and which will be transported to the Gran Sasso Laboratory tunnel as soon as the 
technical teals, are complete. In 1999 a proposal was submitted for the ICANOE 



experiment ,E2!J which combined the ICARUS technology with a fine-grained iron 
calorimeter developed by the NOE collaboration.lli2l 

Because of the high resolution on the event kinematic variables, ICANOE could 
search for Vr appearance in the NCGS beam using kinematic selection criteria, 
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Fig. 26. ICANOE expected visible energy distribution for events with a primary electron, in the 
casp-oF i/ii — > Ut oscillations with maximum mixing angle and Am-^ = 3.5 X 10~^ eV^ (from 
ref.E^). Both the background from the prompt fe and in the beam (hatched histogram) and 
the events from t — » e decays (dashed line) have been considered. Errors are only statistical. 



a strategy which has been successfuhy apphed by the short-basehne experiment 
NOMAD 



112 



For ICANOE, as for NOMAD, the r ^ e decay mode represents a 
golden channel. After four years of running with intensity 4.5 x 10^^ pot/year, the 
expected number of j^t-CC events followed by a t ^ eviy decay would be about 110 
for Am^ = 3.5 x 10"'^ eV^, while the total expected number of background events 
from prompt and is about 470. This means that an excess of events should 
be measured in the visible energy distribution even before applying any kinematic 
cut (see Fig.p6|) . 

The excellent electron identification in ICANOE would also allow a search for 
i^e appearance from — s- oscillations. Selection criteria based on kinematic 
variables could be used to discriminate between events from z^^ — > v,, and events 
from T — > e decays. Assuming Amis = 3.5 x 10~^ eV^ and 623 = 45*^, the statistical 
significance for a — > z/g signal (after four years of running w ith intensity 4.5 x 
10"'^^ pot/year) varies between 0.8 a and 6.8 a for 3*^ < 6*13 < 9*^. 



In summer 2000 the ICANOE collaboration ceased to exist as such, but it is likely 
that as soon as the successful operation of the T600 module is demonstrated, ^11 



updated version of the proposal will be submitted to the relevant organisations .Eli 



Until now we have discussed the experimental hints for neutrino oscillations 
coming from the study of naturally produced neutrinos and we have described the 
experimental program which, in the near future, should clarify many of the open 
issues concerning the solar and the atmospheric neutrino anomalies. There is a 
third hint for neutrino oscillations, which has the peculiarity of being the only one 
coming from an accelerator experiment and the only one to represent evidence for 
neutrino oscillations in the appearance mode. This is the so-called LSND effect, 
which we shall discuss in the next section. 
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7. The LSND EfTect 

The LSND experimenJlii at the Los Alamos Meson Physics Facihty (LAMPF), 
New Mexico, USA, was originally designed to search for i/^ oscillations using 

i/^ from ^+ decay at rest (DAR). In 1995 LSND reported an excess of electron- 
like events above background and gave an interpretation of that result in terms 



of neutrino oscillations.lii3 The effect wa s la ter confirmed by an improved analysis 
with larger statistics, both in DAR z?^'ali3 and in the i^fj, sample from tt"*" decay 



atist 



in flight (DIF)Jiij In the followin g we shall discuss the results reported in a recent 



paper of the LSND collaboration,Ei3 where DAR and DIF events are treated in the 
same way, in a unified analysis of the two samples. 

The neutrino beam used by LSND was obtained from the decay of the secondary 
particles, mostly 7r+ with a small fraction of tt" [N^- /N^+ ^ 1/8), produced in the 
interaction of an intense (~ 1 mA) 798 MeV proton beam with the primary target. 

The decay chains leading to the production of neutrinos and anti-neutrinos are: 

7r+ -^ii+v^ (50) 
^ e^VeV^i (51) 

and the charge-symmetric sequence: 

vr^ -^ji^v^ (52) 
'-^e~VeVf^. (53) 

The fractions of 7r+ DAR and DIF are 96.6% and 3.4% of the total respectively. 
While only ~ 0.05% of the positive muons produced in (50) decay in flight, most of 
them stop before decaying and then give a standard Michel spectrum for v,, and 
{Emax{v^j) = — 52.8 MeV). About 95% of the 7r~'s are stopped in the target 

and immediately captured, while the remaining 5% decay in flight; moreover, about 
88% of the muons from the 7r~ DIF are captured from atomic orbit. In conclusion 
a suppression of the order of - 1/8 x 0.05 x (1 - 0.88) = 7.5 x 10"'' is expected 
for the relative flux of prompt i/g's, which constitute a background for the 
oscillation search. The v^, flux from /i+ and 7r+ DIF, representing a background for 
the v^j_ — > Vf, oscillation search, is suppressed by the long muon lifetime and by the 
small 7r+ e+i^e branching ratio (1.23 x 10^**). 

In 1996 the production target, which was made of water during the running 
period 1993-1995, was upgraded: two upstream targets, used to provide beam to 
other experiments, were removed and the water production target was replaced by 
tungsten. This resulted in a reduction of the DAR flux per pot and of the DIF flux 
per pot by 27% and 34% respectively. 

The estimated systematic uncertainty on the beam predictions is 7% for /x"*" 
fluxes and 15% for ■k'^ DIF and [i' DAR fluxes. 

The LSND detector was an approximately cylindrical tank (8.3 m long, 5.7 m 
diameter) , located 30 m downstream from the neutrino source (sensitivity to Am^ ~ 
0.1—1 eV^) and filled with 167 tons of mineral oil doped with scintillator (0.031 g/liter 
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of butyl-PBD). Because of the low scintillator concentration, both scintillation and 
Cherenkov light could be detected by the 1220 photo-multiplier installed on the 
inside surface of the tank (25% coverage). The Cherenkov cone and the time distri- 
bution of the light were used to tag electron events and to measure the event vertex 
and the direction. For any PMT activity above threshold, the pulse height, a 
digitised time and the event position were recorded, so that energy measurements 
and space and time correlation were possible. On all sides except the bottom the 
detector was surrounded by an active shield (scintillator plus lead), providing a veto 
system against charged particles (10^'' inefficiency). Triggers were not required to 
be in-spill, but the beam state was recorded by the data acquisition and a beam- 
on/beam-off comparison was used to perform a statistical subtraction of the residual 
cosmic ray background. 

The identification of events in LSND is based on the detection of quasi-elastic 
scattering: 

Ue+p^ e'^ +n (54) 

followed by the neutron capture process: 

n+p^ d + ^{2.2MeV) (55) 

Therefore the Pe signature in LSND is a delayed coincidence between an electron 
signal and a spatially correlated 2.2 MeV photon. 

Electron neutrinos are identified in LSND via the inclusive charged-current re- 
action: 

UeC e'X . (56) 

Although the primary electron selection criteria have been optimize for the DAR 

analysis, the same cuts are used to select also electrons for the DIF analysis. The 
electron energy range defined for the v,, search is 20 < i^e < 60 MeV, while 

events with 60 < < 200 MeV are used for the analysis. The cut at 

20 MeV rejects the background events from the /?-decay of ^^B produced in the 
capture of cosmic n~ on ^^C, while events above 200 MeV are removed to suppress 
the beam-related background from tt"*" e+i'e, which at those energies becomes 
larger that any possible oscillation signal. Cosmic ray events are further suppressed 
by using the veto information and by removing electron events with additional 
activity within 8 /xs before the electron signal or 12 /xs after it. Discrimination 
between correlated 2.2 MeV photons from neutron capture and accidental 7's from 
radioactivity is achieved by means of a likelihood function, R^^ which depends on 
the number of PMT hits, the distance between the reconstructed 7 and e"*" positions 
and the time interval between the 7 and the e+ events. The discrimination power of 
has been checked by using Vf,C — > e~Ng,s, exclusive interactions, where no recoil 
neutron is present and which should therefore exhibit an R^ distribution consistent 
with that expected for purely accidental 7's. 
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The two largest backgrounds to the i?^ signal are VeCG interactions from 

prompt Vf, in the beam, which are identical to the signal events, or v^{v^)CC inter- 
actions with a correlated neutron signal and the muon misidentified as a positron. 
Both backgrounds have been estimated by using the experiment Monte Carlo sim- 
ulation and, whenever possible, by means of control samples from real data. 

For the entire period 1993-1998 of data-taking, LSND observed a beam on-ofF 
excess of 117.9±22.4 events in the DAR channel. The corresponding expected num- 
ber of background events from /i^ DAR followed by v^p —* e+n scattering and from 
7r~ DIF followed by i/^p fX^n scattering are 19.5 ±3.9 and 10.5 ±4.6 respectively, 
leading to a total excess of 87.9±22.4±6.0 events above background expectationst 
The number of expected Ve events from — > Ve oscillations with 100% probability 
is 33000 ±3300, the error being a combination of the systematic uncertainties on the 
neutrino flux (7%) and on the e+ efflciency (7%). The result then corresponds to an 
oscillation probability of (0.264 ± 0.067(siai.) ± 0.045(sys.))%, to be compared 
with the previously published result of (0.31 ± 0.12(stat.) ± 0.05(s?;s.))%, 1113 which 
was based on the 1993-1995 data sample. 

The new DIF v^. oscillation analysis gives a beam on-off excess of \^.l:kY1.2 

events. The estimated total background from the /i^ — s- e^i'ei^^, tt"*" — s- and 
ve ve processes is of 6. 6± 1.7 Vf^ events, resulting in a total excess of8.1±12.2±1.7 
events above background expectations. The number of events expected for v^^ — s- v^, 
oscillation with 100 % pr obability is 7800, thus giving an oscillation probability of 
(0.10±0.16±0.04)%.ti3 This has to be compared with the result obtaine d fro m the 



1993-1995 data analysis, which gave a probability of (0.26 ± 0.10 ±0.05)%H3 Using 
the result from the DAR analysis and assuming that CP is conserved, the expected 
oscillatio n pr obahility in the DIF sample is ~ 0.26% at high Am^ and ~ 0.05% at 
low Am2.lii 

A likelihood fit of the beam-on events which pass their oscillation cuts is used 
to extract the LSND allowed region in the sin^ 20, Am^ plane. The fit is performed 
over the entire electron range 20 — 200 MeV, so that both DAR and DIF events 
are considered at the same time. The probability density functions which define 
the likelihood function £ used in the fit depend on the electron energy, the electron 
longitudinal position along the tank, the cosine of the reconstructed angle between 
the neutrino and the electron direction, the likelihood ratio R-^ for correlated 7 
selection. To account for the limited knowledge of the background, its contribution 
to £ is allowed to vary with gaussian profile around the expected central value." 

The result obtained is shown in Fig.|2^, where the LSND allowed region is com- 
pared with limits obtained by other experiments sen sitive to v, , --> v,, {vn — > Vg) 



oscillations in different regions of the parameter space. E^illtlilHiij'lirj In particular it 



is interesting to briefly discuss the null result obtained by the KARMEN2 collabo- 
ration, which has conducted a search for v^ oscillations in approximately the 
same region of the parameter space as that investigated by LSND. 



°beam-related background are locked to vary together, while beam-unrelated backgrounds can 
fluctuate independently. 
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The KARMEN2 experiment was performed at the ISIS neutron spallation facil- 
ity at the Rutherford Appleton Laboratory, Chilton, Didcot, United Kingdom. As 
at LAMPF, neutrinos are obtained from the decay at rest or in flight of secondary 
particles produced in the interactions of 800 MeV primary protons with a massive 
beam stop target. Also in this case the I'^'s from n'^ decays are monoenergetic 
{E^^ = 29.8 MeV), while both i/e's and i?^'s have continuous energy spectra ex- 
tending up to 52.8 MeV. One difference between the two beams is that the one at 
ISIS has a well defined time structure: two very short 100 ns) proton pulses, 
separated by 225 ns, are produced with a repetition rate of 50 Hz. Because of the 
different lifetimes of their parent particles (ttt — 26 ns, = 2.2 /is), the ^'^-burst, 
which occurs essentially in coincidence with the proton pulses, can be clearly sep- 
arated from the z/g- and //^-induced events, which are instead characterised by a 
much longer decay constant. Due to the accelerator duty cycle of 10~^ (to be com- 
pared with the 6 x 10^^ duty factor for LSND), the cosmic ray background is largely 
suppressed with respect to the beam-related events. The proton beam intensity at 
ISIS is 0.2 mA, that is to say a factor of five lower than at LAMPF. 

The KARMEN detector,E^ located 17.5 m downstream the neutrino source 
(about half the LSND baseline), consisted of a rectangular tank filled with 56 t 
liquid scintillator. The tank was segmented into 512 optically independent modules, 
wrapped with Gd203 coated paper to increase the efficiency for the neutron capture 
process. The event position was obtained from the time of the PMT signal at the two 
ends of the hit module. Shielding against the beam-related neutron background and 
the cosmic ray background was provided by a 7 kton steel blockhouse combined with 
a system of two layers of veto counters. Moreover an additional veto counter, added 
in 1996 for the KARMEN2 upgrade, was used to largely suppress the background 
from energetic neutrons produced in I'pCC DIS interactions. 

The De detection in KARMEN2 was based on the delayed coincidence between 
the positron signal from reaction ( p4|) and either the 2.2 MeV photon of reaction 
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( psf ) or the 8 MeV integrated photon signal from Gd de-excitation following neutron 
capture on gadolinium. No particle identification was performed and the tagging 
of the positron event, which was expected within 2.2 /is after the beam-on-target 
signal, relied on the beam time structure alone. Both time and spatial correlation 
between the positron and the neutron capture events were required. An analysis 
of the data collected between February 1997 and March 2000 shows that 11 events 
have been, identified as De interactions, with an expected background^ of 12.3 ± 0.6 
events. E£j Therefore the KARMEN2 result is compatible with the no-oscillation 
hypothesis and the corresponding exclusion plot at 90% C.L. is shown in Fig.^. 

Combining the LSND result with those obtained by all the other experiments 
which, having studied the same transition, have obtained null results, only a small 
portion of the LSND solutions is still allowed to explain the experiment anomaly 
in terms of neutrino oscillation. In particular, the skeptical reader may find the 
LSND and KARMEN2 results quite contradictory. However it has to be stressed 
that the two experiments have sensitivities which peak at slightly different values 
of Am^. Moreover, before any definitive conclusion can be drawn from the various 
results, a combined analysis of all experimental data should be performed. This has 
been tried for LSND and KARMEN, but only on subset of the collected dataJHj 
In conclusion it is fair to say that no definitive answer has yet been given to the 
LSND puzzle and that more expe rime ntal data is needed. 

The MiniBOONE experiment ,E£l which is now being built at Fermilab, wil 
address the LSND effect by searching for z^^ Vf, oscillations in the same region 
of parameter space. An almost pure beam (Vf, contamination < 0.3%), with 
energy in the range 0.5 — 1 GeV, will be obtained from the decay of secondary 
particles produced in the interactions of 8 GeV protons with a beryllium target 
positioned in a magnetic horn. The protons will be delivered by the Fermilab 
Booster, with an expected intensity of 5 x 10^" pot per year. The decay length for 
the positively charged secondary mesons, which are focused by the single horn, can 
be changed from 50 m to 25 m by inserting an intermediate absorber inside the 
decay tunnel, thus allowing a systematic study of the prompt v^. contamination of 
the beam. The MiniBOONE detector will be located about 500 m from the neutrino 
production point, so that the E j L ratio will be approximately the same as LSND's. 
MiniBOONE will consists of a spherical tank (12.2 m diameter) filled with 800 tons 
of pure mineral oil. The Cherenkov light produced in neutrino interactions in the 
inner part of the detector (445 ton fiducial volume), as well as the modest amount 
of scintillation light from intrinsic impurities in the oil, will be detected by 1280 
PMTs (10% coverage) installed on a support structure in the tank. This structure 
will also serve as an optical insulator for the inner part of the detector, while the 
outer 35 cm thick oil layer, read out by 280 additional PMTs, can be used as a veto 
system. 

The main background components will be due to the v^. contamination in the 



P Background events come from cosmic muon interactions, y]rC — > e ^g.a. sequences, prompt 
v„ contamination in the beam and neutrino-induced accidental coincidences. 
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beam and to the misidentification of muons and vr*^ as electrons.^ The estimated 
systematic uncertainties on the i/^. flux is about 5% for the neutrino component 
originating from the 7r+ fj,^ —^ chain and 10% for v^s coming from iiTeS 
decays of secondary K'^ and K'j^. Particle identification will be performed using 
the sharpness of the Cherenkov rings and the ratio between the amount of prompt 
(Cherenkov) and late (scintillation) light, which is different for electrons and muons. 
The misidentification uncertainty is expected to be less than 5% for muons and 
about 5% for neutral pions: in the case of muons, the the systematic error can be 
estimated using a sample of muon decays from real data, while for tt'^'s it will be 
obtained from a Monte Carlo simulation, which can in turn be constrained using 
reconstructed 7r°'s from real data. 
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Fig. 28. Final 90% C.L. expected sensitivity for MiniBOONE (figure from ref.E£]). Different beam 
configurations and cut efficiencies have been considered. 



If oscillations occur with parameters compatible with LSND solutions, after one 
year of running MiniBOONE should observe an excess of ^ 1500 Vf, events above 
background expectations, with a significance of ~ 8 — lOcr. On the other hand, 
in case of null result, after two years of running MiniBOONE should be able to 
completely exclude the entire LSND allowed region at 90% C.L. (see Fig.|2^). 

8. Neutrino Masses and Global Analysis of Neutrino Oscillation Data 

If we believe that all experimental hints on neutrino oscillations are correct and 
we accept that they are all manifestations of flavour mixing in the lepton sector, 
we are forced to introduce at least another neutrino species in our framework. In 
fact, the Am? regimes suggested by the the solar neutrino deficit, the atmospheric 
neutrino anomaly and the LSND effect are all well disconnected one from the other: 

Am| « Aml^,^ « AtjiIsnd ■ (57) 

Given the higher neutrino energy with respect to LSND, neutrons will not constitute a significant 
background for MiniBOONE 
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On the other hand, only two independent values of Am^ are possible for three 
different neutrino masses [Arn^i = Aml2 + Amfi). Therefore, recalling the LEP 
result on the number of active neutrinos,El in order to accommodate all experimental 
results, an additional sterile neutrino would be required. 

While both the solar and the atmospheric neutrino deficit have been confirmed 
by more than one experiment, using different experimental techniques, no other 
experiment has confirmed LSND observations. For this reason, given also the enor- 
mous impact that accepting the LSND result would have on our understanding of 
sub-nuclear phenomena, with the introduction of a light sterile neutrino, there is 
a tendency in the scientific community to adopt a "conservative" attitude and to 
exclude LSND from most phenomenogical models and fits to neutrino oscillation 
parameters, until more solid experimental data become available. In the following 
we shall take the same approach and discuss some of the results obtained from three- 
neutrino analyses of the neutrino oscillation data, while only briefly commenting on 
the possibility to have a fourth neutrino. 

Neglecting possible CP violation phases, which are anyhow not accessible to 



present experiments, the lepton mixing matrix (13) depends only on five parame- 
ters, two Am^ and three mixing angles, four of which can be identified with the 
parameters measured by solar and atmospheric experiments: 

Atoq = Am22 
Am^,„ = Amj^ (58) 

Oatm = 6*23 

while a limit on the fifth, ^13, is set by the CHOOZ reactor experiment. 

The hierarchy Auiq ^ Am^(,„ implies that one of the three mass eigenstates, 
say 1^3, is separated by a larger gap |Arn|;^| w |Am§2| ~ Am^^^ ~ 10~^ eV^ from 
the other two, vi and V2, which are instead quasi-degenerate, being separated by 
jArrijil ~ Atoq < 10^"' eV^. However, since oscillations determine only mass 
splittings and not the absolute value of the masses, it is still possible to shift the 
absolute mass scale without affecting the oscillation phenomenology. Although dif- 
ferent mass orderings are possible, here we consider the normal hierarchical scheme, 
corresponding to: 

mi < 7712 <C rn^ (59) 

and thus: 

Arn2i < Am32 w Arng^ . (60) 

It can be shown that, given the strong hierarchical constraint Atoq <C Am^^^, 
the solar neutrino data can be described in terms of only three parameters, Am^i, 
9i2 and ^13, while the atmospheric neutrino data can be analysed in terms of Am§2, 
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623 and 6*13. For the CHOOZ data, the survival probabihty depends only on 6*13 
and Am§2- Thus ^13 is the only parameter common to the three data sets. In 
particular, in the limit ^13 = atmospheric and solar neutrino oscillation can be 
described by two independent two-flavour oscillation analyses. 

SevemLanalyses of the available data have been performed in the three-neutrino 
schemeJ2d'E£3 in the following we shall refer to the results presented in ref.E2l. In 
particular it will be interesting to see how the inclusion of different data in the fit 
affects the constraint which can be extracted on the mixing angle ^13. 

The three-dimensional analysis of ref.L3 of the solar neutrino data has shown 
that solutions with small values of the mixing angle ^13 are preferred, the most 
favourable scenario being that of two-flavour mixing (^13 = 0). The upper limit 
which can be extracted on the mixing angle based on solar neutrino data only is 
tan2 6'i3 < 2.4 (6*13 < 57") at 90% C.L.. 

If only atmospheric neutrinos are considered, the analysis favours the i'^ Vt 
oscillation hypothesis. Also in this case the best fit corresponds to a small value of 
the mixing angle 6*13 (Q\3 — 9"). However oscillations between muon and electron 
neutrinos are still allowed as a sub-dominant mode, with the following 90% C.L. 



allowed ranges for the relevant oscillation parameters: 

1.6 X 10"^ eV^ <Am^2 < 6 x 10"^ eV^ (61) 

0.43 <tan2 6*23 < 4.2 (62) 

tan^ 6*13 < 0.34 (63) 



If the CHOOZ result is included in the atmospheric neutrino analysis, because 
the full allowed region of Am|2 lies within the parameter domain accessible to 
CHOOZ, a tighter constraint can be obtained on the mixing angle ^13 and the 
corresponding 90% C.L. limit becomes even more stringent: tan^ ^13 < 0.043 (6*13 < 
12°). 

Finally, a combined analysis of solar, atmospheric and reactor data gives, at 
90% C.L., the following allowed ranges for Am§2, tan^ ^23 and tan^0i3, where the 
limits on tan^ 6*13 are given for the SMA solution of the solar neutrino deficit as well 



as for the LMA or unconstrained case:'' 

1.4 X 10"^ eV^ <Am^2 < 6 x 10"^ eV^ (64) 

0.39 <tan2 6*23 < 3.1 (65) 

tan^ 6*13 < 0.055 (unconstrained or LMA) (66) 

tan^ 6*13 < 0.075 (SMA) (67) 



The allowed regions for Ato2i and 6*12 obtained from this global analysis are 
shown in Fig.^, both for the unconstrained and constrained case. It is interesting 
to notice that, in the unconstrained case, the best-fit solution lies in the LMA 
region. 

' in the unconstrained case no assumption is made about the solution to the solar neutrino deficit. 



58 An Experimentalist's View of 




IS* 10'^ ler^ 10'^ is^ 10'^ i lo. lo"* lo'^ iff^ lo'^ ler^ ler^ i lo. loo. 

Fig. 29. Allowed regions in the plane {tan^ 612, Amju^rias obtained from the global three-neutrino 
analysis of solar, atmospheric and reactor data of refH, for both the unconstrained (left) and the 
constrained fit (right). The star denotes the global best fit solution for the unconstrained fit and 
local best fit points for the constrained fit. 



In summary, the analysis in ref.0 shows that, although all solutions to the 
solar neutrino deficit are still allowed, the LMA solution appears to be favoured. 
Moreover both solar and atmospheric neutrino data seem to prefer small values of 
the mixing angle ^13, the constraint becoming more stringent with the inclusion of 
the CHOOZ result. Limits have been extracted also on the other relevant oscillation 
parameters. 

Ignoring CP violation, in the limit Ues ^ (that is to say S13 ~ and cia ~ 1) 
the MNS matrix (O) can be written as: 



C12 S12 
U ~ I -S12C23 C12C23 
S12S23 — C12S23 



(68) 



One appealing realization of this ma trix would be that corresponding to bi- 



maximal mixing (see, for example, ref.ti3), which postulates that both atmospheric 



and solar neutrino anomalies can be interpreted as neutrino oscillations with max- 
imal mixing (that it to say both 623 and 612 equals 7r/4). In this case: 



(69) 
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However it has to be noticed that this form of the mixing matrix seems to be 
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incompatible with the LMA solution of the solar neutrino deficit .li^ 

If one tries to define a coherent picture which accommodates also the LSND 
result, a fourth neutrino has to be introduced. Essentially only two mass schemes 
can used to describe the mass spectrum, the so called 3 + 1 and 2 + 2 sc hemes. The 



3+1 scheme, which seems to be disfavoured by experimental data,E3'E£j corresponds 
to having three relatively close neutrino masses, whose spacing is determined by 
the two Am^ from the solar and atmospheric neutrino analyses, separated by about 

1 eV^ from an isolated level, which is responsible for the LSND anomaly. In the 

2 + 2 scheme, instead, two doublets of mass eigenstates, spaced by Atoq and Am^j^ 
respectively, are separated by the LSND gap. Several analyses of the solar and 
atmospheric neutrino data in the 2 + 2 four-neutrino framework have been recently 
presented, showing that, although Super-Kamiokande data disfavour both the pure 



and the pure — > chauni 



with sizable oscillation amplitude 



9. Conclusions 



133, L34 



fi^t is not possible to exclude their occurrence 



In this paper we have reviewed the current status of the search for neutrino 
oscillations. After an introduction of the formalism of neutrino oscillation phe- 
nomenology, we have discussed the various experimental results which, over the 
past three decades, more or less convincingly, have provided hints for the existence 
of flavour mixing in the neutrino sector: the solar neutrino deficit, the atmospheric 
neutrino anomaly and the LSND effect. If any of these indications will be conflrmed 
by future experiments, we would be forced to extend our description of sub-nuclear 
phenomena beyond the Standard Model. 

Having described the different experimental techniques used to explore different 
domains of the oscillation parameter space, and having discussed the sources of 
systematic error intrinsic to each method, we hope we have given most of the in- 
formation necessary to discern the robustness of each method and result. However, 
whatever the preference of the reader may be, it is undeniable that more exper- 
imental data is urgently needed, before a clearer picture of neutrino oscillation 
phenomenology can be drawn. 

The near and medium term experimental programme, which is supposed to 
answer many of the as yet unresolved puzzles, has been described in the previous 
sections. However it is possible that the experiments already approved will not be 
able to provide a definitive solution to all our pending questions. In particular, none 
of the experiments scheduled to start running over the next few years will have the 
sensitivity to address CP violation in the neutrino sector and, except if the actual 
value of 6*13 is very close to the present CHOOZ limit, none of them will be ableJx 
measure that mixing angle. Very ambitious pr ojec ts, such as neutrino factories. 



or super-intense neutrino beams [super-beams) if ever approved, might hold the 
key to the ultimate solution to the fascinating neutrino mystery. 

To conclude, whatever the outcome of future experiments, there is no doubt this 
is an extremely exciting time for neutrino physics. 
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